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*^ In these lectures we discuss the application of discrete light cone quantiza- 

tion (DLCQ) to supersymmetric field theories. We will see that it is possible to 
formulate DLCQ so that supersymmetry is exactly preserved in the discrete ap- 
proximation. We call this formulation of DLCQ, SDLCQ and it combines the 
power of DLCQ with all of the beauty of supersymmetry. In these lecture we will 
review the application of SDLCQ to several interesting supersymmetric theories. 
We will discuss two dimensional theories with (1,1), (2,2) and (8,8) supersymme- 
try, zero modes, vacuum degeneracy, massless states, mass gaps, theories in higher 
dimensions, and the Maldacena conjecture among other subjects. 
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Introduction. 

In the last decade there have been significant improvements in our understanding of 
gauge theories and important breakthroughs in the nonperturbative description of su- 



persymmetric gauge theories |74], [75|]. In the last few years various relations between 
string theory, brane theory and gauge fields |3^, 0] have also emerged. While these devel- 
opments give us some insight into strongly coupled gauge theories |7Sj , they do not offer a 
direct method for non-perturbative calculations. In these lectures we discuss some recent 
developments in light cone quantization approaches to non-perturbative problems. We 
will see that these methods have the potential to expand our understanding of strongly 
coupled gauge theories in directions not previously available. 

The original idea was formulated half of a century ago |32|, but apart from several 



technical clarifications f63[ it remained mostly undeveloped. The first change came in 



the mid 80-th when the Discrete Light Cone Quantization (DLCQ) was suggested as a 
practical way for calculating the masses and wavefunctions of hadrons [ 68| . Although 



the direct application of the method to realistic problems meets some difficulties (for 
[25|j), DLCQ has been successful in studying various two dimensional models. 



review see 



Given the importance of supersymmetric theories, it is not surprising that light cone 



quantization was ultimately applied to such models |57], |20], p23f| . Even in this early 
work the mass spectrum was shown to be supersymmetric in the continuum and a great 
deal of information about the properties of bound states in supersymmetric theories was 
extracted. However the straightforward application of DLCQ to the supersymmetric 
systems had one disadvantage: the supersymmetry was lost in the discrete formulation. 



The way to solve this problem was suggested in |64| , where an alternative formulation of 
DLCQ was introduced. Namely it was noted that since the supercharge is the "square 
root" of Hamiltonian one can define a new DLCQ procedure based on the supercharge. 
We will study this formulation (called SDLCQ) in these lectures. 

These lectures have the following organization. In section [l] we introduce the basic 
concepts of DLCQ and SDLCQ. We also define the systems to be studied in the remaining 
lectures. We will concentrate our attention on two dimensional models with adjoint 
matter, several examples of such systems can be constructed from supersymmetric Yang- 
Mills theories in higher dimensions using dimensional reduction. In section we address 
the problem of the DLCQ vacuum. In the continuum theory the light cone vacuum is 
very simple: it coincides with the usual Fock vacuum. This property is related to the 
decoupling between positive and negative frequency modes on the light cone and does 
not occur for equal time quantization. In DLCQ however one encounters the problem of 
zero modes which complicates the structure of vacuum and allows us to reproduce the 
correct vacuum degeneracy in certain theories. We continue to analyze zero modes in 
the section [3] and they are shown to play an important role in explaining the difference 
between DLCQ and SDLCQ regularization procedures. 

Section [| is devoted to the study of massless states. Our numerical analysis |M|, 
H [11, |J shows an important property of the mass spectrum for some supersymmetric 



theories. We find that unlike the QCD-like models [EH], such systems appear to have a 



lot of massless bound states and in fact the supersymmetric SU(oo) gauge theory seems 
to have an infinite number of such states in the continuum limit. Since the states with 
zero mass dominate the partition function for low enough temperatures they deserve to 
be studied very carefully and in section |] we analyze the structure of such states. 

As we already mentioned in the beginning of this introduction, the relation between 
string theory and gauge fields has attracted a lot of attention in recent years. In particu- 
lar it was conjectured |6l| that one can extract some information about strongly coupled 



gauge theory from supergravity calculations. The problem however is that in the rele- 
vant regime usual field theoretic methods do not work, so it is hard to really test the 
conjecture. For two dimensional systems, however DLCQ gives solutions of the bound 
state problem which are valid beyond perturbation theory, so the results can be used to 
test the conjecture. We report the results of this first test in section [5]. For realistic sys- 
tems with eight supersymmetries we still don't have enough computer power to compare 
the results with the supergravity predictions. The general techniques described in this 
section can also be used to calculate other correlation functions in the nonperturbative 
regime. 

Finally in section |6] we make a first attempt to move beyond two dimensions. We 
present the general ideas for formulating SDLCQ in more than two dimensions. As an 
example we present the numerical results of SYM for the simplest case when only one 
transverse momentum mode is introduced. 

1 Supersymmetric Yang— Mills Theory in the Light- 
Cone Gauge. 

1.1 DLCQ and Its Supersymmetric Version. 

In this work we will study the bound state problem for various supersymmetric matrix 
models in two dimensions. The examples of such models may be constructed by di- 
mensional reduction of supersymmetric Yang-Mills theory in higher dimensions. In this 
subsection we will consider such reduction for three dimensional SYM. Before we begin 
the detailed analysis of the bound state problem for the specific systems it is worthwhile 
to summarize some basic ideas of Discrete Light Cone Quantization ( for a complete 
review see |[25|| ) . 

Let us consider a general relativistic system in two dimensions. Usual canonical 
quantization of such a system means that one imposes certain commutation relations 
between coordinates and momenta at equal time. However as was pointed out by Dirac 
long ago [B3] this is not the only possibility. Another scheme of quantization treats the 
light like coordinate x + = 75(2° + x 1 ) as a new time and then the system is quantized 
canonically. This scheme (called light cone quantization) has both positive and negative 
sides. The main disadvantage of light cone quantization is the presence of constraints, 



even for systems as simple as free bosonic field. From the action 

S = f d 2 xd + (j)d-(j) (1.1) 

one can derive the constraint relating coordinate and momentum: 

7T = <9_0. (1.2) 

For more complicated systems the constraints are also present and in general they are 
hard to resolve. 

The main advantage of the light cone is the decoupling between positive and negative 
momentum modes. This property is crucial for DLCQ. In the Discrete Light Cone Quan- 
tization one considers the theory on the finite circle along the x~ axis: — L < x~ < L. 
Then all the momenta become quantized and the integer number measuring the total 
momentum in terms of "elementary momentum" is called the harmonic resolution K. 
Due to the decoupling property one may work only in the sector with positive momenta 
where there are a finite number of states for any finite value of resolution. Of course the 
full quantum field theory in the continuum corresponds to the limit L — > oo and , in this 
limit the elementary bit of momentum goes to zero, as the harmonic resolution goes to 
infinity and the infinite number of degrees of freedom are restored. However it is believed 
that the " quantum mechanical" approximation is suitable for describing the lowest states 
in the spectrum. Note that the problem of constraints in DLCQ is a quantum mechanical 
one and thus it is easier to solve. Usually this problem can be reformulated in terms of 
zero modes and the solution can be found for any value of the resolution. 

DLCQ is mainly used in order to solve the bound state problem. Let us formulate this 
problem for general two dimensional theory. The theory in the continuum has the full 
Poincare symmetry, thus the states are naturally labeled by the eigenvalues of Casimir 
operators of the Poincare algebra. One such Casimir is the mass operator: M 2 = P^P^ 



A'' 



Another Casimir is related to the spin of the particle and we will not use it. After 
compactifying the x~ direction one looses Lorenz symmetry, but not the translational 
invariance in x + and x~ directions. Thus P + and P~ are still conserved charges, but 
now the mass operator is not the only Casimir of the symmetry group: the states are 
characterized by both P + and P~ . However if we consider DLCQ as an approximation 
to the continuum theory we anticipate that in the limit of infinite harmonic resolution 
(or L — > oo) the Poincare symmetry is restored and the mass will be the only quantity 
having invariant meaning. Thus the aim would be to study the value of M 2 as function 
of K and to extrapolate the results to the K = oo. 

The usual way to define M 2 in DLCQ is based on separate calculation of P + and P~ 
in matrix form and then bringing them together: 

M 2 = 2P + P-. (1.3) 

Usually one works in the sector with fixed P + , and the calculation of light cone Hamilto- 
nian P~ is the nontrivial problem. An important simplifications occur for supersymmetric 



theories I64 



Supersymmetry is the only nontrivial extension of Poincare algebra compatible with 



the existence of the S matrix ||80| . Namely in addition to usual bosonic generators of 
symmetries, fermionic ones are allowed and the full (super)algebra in two dimensions 
reads: 

{QlQi} = 25 IJ 1 ^P, + e aP Z IJ , (1.4) 

[P»,Pu] = 0, [P^Qi] =0. (1.5) 

In this expression e is an antisymmetric 2x2 matrix, e 12 = 1 and Z IJ is the set of 
c-numbers called the central charges. In these lectures we will put them equal to zero. It 
is convenient to choose two dimensional gamma matrices in the form: 7 = a 2 , 7 1 = ia 1 , 
then one can rewrite ( |1.4p in terms of light cone components: 

{QhQ + j} = 2V2S IJ P+, (1.6) 

{QJ,Qj} = 2V28 IJ P~, (1.7) 

{Qf,Qj} = 2Z U . (1.8) 

As we mentioned before, in DLCQ diagonalization of P + is trivial and the construction 
of Hamiltonian is the main problem. The last set of equations suggests an alternative 
way of dealing with this problem: one can first construct the matrix representation 
for the supercharge Q~ and then just square it. This version of DLCQ first suggested 
in [64] appeared to be very fruitful. First of all it preserves supersymmetry at finite 
resolution, while the conventional DLCQ applied to supersymmetric theories doesn't (we 
will consider the relation between these two approaches in section 0). The supersymmetric 
version of DLCQ (SDLCQ) also provides better numerical convergence. 

To summarize, in this subsection we defined two procedures for studying the bound 
state spectrum: DLCQ and SDLCQ. To implement the first one we construct the light 
cone Hamiltonian and diagonalize it, while in the second approach one constructions the 
supercharge and from it the Hamiltonian. Of course the SDLCQ method is appropriate 
only for the theories with supersymmetries, although it can be modified to study models 
with soft supersymmetry breaking (see section |6]). 

1.2 Reduction from Three Dimensions. 

Let us start by the defining a simple supersymmetric system in two dimensions. It can be 
constructed by dimensional reduction of SYM from three dimensions to two dimension. 
The more general case can be found in the next subsection. 

Our starting point is the action for SYM in 2 + 1 dimensions: 

S = f d 3 xti (-jF AB F AB + i^ A D A q) . (1.9) 

The system consists of gauge field A A and two-component Majorana fermion \1/, both 
transforming according to adjoint representation of gauge group. We assume that this 



group is either U(N) or SU(N) and thus matrices Af, and *&ij are hermitian. Studying 
dimensional reduction of SYMp we introduce the following conventions for the indices: 
the capital latin letters correspond to D dimensional spacetime, greek indices label two 
dimensional coordinates and the lower case letters are used as matrix indices. According 
to this conventions the indexes in ( |1.9| ) go from zero to two, the field strength Fab and 
covariant derivative Da are defined in the usual way: 

Fab = d A A B - d B A A + ig[A A , A B ], 

D A V = d A * + ig[A A ,% (1.10) 

Dimensional reduction to 1 + 1 means that we require all fields to be independent 
on coordinate x 2 , in other words we place the system on the cylinder with radius L± 
along the x 2 axis and consider only zero modes of the fields. The possible improvement 
of this approximation will be suggested in section || here we consider this reduction 
as a formal way of getting two dimensional matrix model. In the reduced theory it is 
convenient to introduce two dimensional indices and treat A 2 component of gauge field 
as two dimensional scalar <p. The action for the reduced theory has the form: 

S = J d 2 x tr (-1f iiv F> h ' + \d^D^ + M-fDfl - 2ig0j 5 ^j , (1.11) 

We also could choose the special representation of three dimensional gamma matrices: 

7° = a 2 , 7 1 = icr 1 , 7 2 = ia 3 , (1.12) 

then it would be natural to write the spinor \P in terms of its components: 

* = (^X) T - (1-13) 

Taking all these definitions into account one can rewrite the dimensional reduction of 
as: 

S = L ± Jd 2 x(^D II (j)Dy + iV2^D + ^ + iV2xD_x+ 
The covariant derivatives here are taken with respect to the light cone coordinates: 

XU ryj -L 

x= = ^. (1.15) 

Note that by rescaling the fields and coupling constant g we can make the constant L± 
to be equal to one, so below we simply omit this constant. 



The bound state problem for the system (|1.14|) was first studied in [64J. The super- 



symmetric version of the discrete light cone quantization was used in order to find the 



mass spectrum. However the zero modes were neglected by authors of | 6"3}j , so we spend 
some time studying this problem in the next section. As we will see, while zero modes 
are not very important for calculations of massive spectrum, they play crucial role in the 
description of the vacuum of the theory. 

Let us consider ( |1.14| ) as the theory in the continuum. In this case one can choose 
the light cone gauge: 

A + = 0, (1.16) 

then equations of motion for A~ and \ give constraints: 

-d 2 _A-=gJ + , (1.17) 

y/2id. x = 9[M], (1-18) 

J + (x) = ^[0(x),«9_0(x)] - -L{V>(x),V(*)}. (1-19) 

Solving this constraints and substituting the result back into the action one determines 
the Lagrangian as function of physical fields <fi and ip only. Then using the usual Noether 
technique we can construct the conserved charges corresponding to the translational 
invariance: 

P+ = [ dx'tr ((d_</>) 2 + iV2i/)d^) , (1.20) 



P- = | cir -tr^J+l-J+ + ^[0 l V]^^]J. (1-21) 

We can also construct the Noether charges corresponding to the supersymmetry trans- 
formation. However the naive SUSY transformations break the gauge fixing condition 
A + = 0, so they should be accompanied by compensating gauge transformation: 

SA, = l -s^ - D^.j-% (1.22) 

The resulting supercharges are: 

Q+ = 2 [ dx'tr {ipd-fl , (1.23) 

Q~ = -2g f dx'ti \J + jrV\ ■ (i- 24 ) 

Finally we make a short comment on supersymmetry in the pure fermionic system. As 
one can see the expression for Q~ contains the term cubic in fermions, so if we formally 
put (j) = this supercharge will not vanish. One may ask what kind of supersymmetric 



system this supercharge corresponds to. This answer was found by Kutasov |57|] who 



discovered the supersymmetry in the system of adjoint fermions, namely the square of 
supercharge including fermions only gives Hamiltonian: 



im 2 1 g 2 . 1 



p- = / dx-ti ( — — v>— i> - ^J + iv J+ > ( L25 ) 



rn 2 = g 2 N/n. This P corresponds to the system of adjoint fermions in two dimensions 
with the special value of mass m. We will consider this system in details in section |3|. 

1.3 Reduction from Higher Dimensions. 

In this subsection we consider the general reduction of SYM^ to two dimensions. By 
counting the fermionic and bosonic degrees of freedom one can see that the SYM can be 
defined only in limited number of spacetime dimensions, namely D can be equal to 2, 
3, 4, 6 or 10. The last case is the most general one: all other system can be obtained 
by dimensional reduction and appropriate truncation of degrees of freedom. So in this 
subsection we will concentrate on reduction 10 — > 2, and the comments on four and six 
dimensional cases will be made in the end. 

As in the last subsection we start from ten dimensional action: 

S= f d 3 xtr (--F AB F AB + fflj A D A y\ . (1.26) 



According to our general conventions the indexes in ( |1.26|) go from zero to nine, \l/ 



is the ten dimensional Majorana- Weyl spinor. A general spinor in ten dimensions has 
2io/2 _ 22 complex components, if the appropriate basis of gamma matrices is chosen 
then Majorana condition makes all the components real. Since all the matrices in such 
representation are real, the Weyl condition 

rn* = * (i.27) 

is compatible with the reality of \l/ and thus it eliminates half of its components. In the 
special representation of Dirac matrices: 

r° = a 2 ® 1 16 , (1.28) 

r / = ia 1 ®7 7 , 1=1,..., 8; (1.29) 

r 9 = iCT!®7 9 , (1.30) 

the r n = r° • • -T 9 has very simple form: r n = a 3 <g) 1 16 . Then the Majorana spinor of 
positive chirality can be written in terms of 16-component real object ip: 



* = 2 1 / 2 r. (i.3i) 




Let us return to the expressions for T matrices. The ten dimensional Dirac algebra 

{r M , rj = 2g lxv 

8 



is equivalent to the spin(8) algebra for 7 matrices: {7/,7j} = 28 u and the ninth matrix 
can be chosen to be 7 9 = 7 1 . . . 7 s . Note that the 16 dimensional representation of spin(8) 
is the reducible one: it can be decomposed as 8 S + 8 C 

7'=(^*). /=!,... ,8. (1.32) 

The explicit expressions for the /3j satisfying {/3i,(3j} = 28 u can be found in p6| . Such 
choice leads to the convenient form of 7 9 : 

7 9 =f !? "I (1-33) 



-1 



N 



So far we have found nonzero components of the spinor given by ( |1.31|) . However 
as we saw in the last subsection not all such components are physical in the light cone 
gauge, so it is useful to perform the analog of decomposition ( [1.131 ). In ten dimension 
it is related with breaking the sixteen component spinor ijj on the left and right-moving 
components using the projection operators 

P L = \{l-l\ P R = i(l+ 7 9 ). (1.34) 

After introducing the light-cone coordinates x ± = -4= (x° ± x 9 ) the action (|1.26|) can be 
rewritten as 



S££ = J dx + dx~dx x tr (~-F+_ + F +I F.j - - 

+ iV2iPlD + 4> R + iV2i>lD„4> L + 2iiPl 1 I D I ip R ), (1.35) 



-/li 



where the repeated indices /, J are summed over (1, . . . , 8). After applying the light-cone 
gauge A + = one can eliminate nonphysical degrees of freedom using the Euler-Lagrange 
equations for ip L and A~: 

d-^L = -^D^n, (1.36) 

d 2 _A + = d_<9Mj + gJ + (1.37) 

J + = i[Aj, d.Aj] + 2^ T R i) R . (1.38) 

Performing the reduction to two dimensions means that all fields are assumed to be 
independent on the transverse coordinates: di§ = 0. Then as in previous subsection one 
can construct the conserved momenta P in terms of physical degrees of freedom: 

P + = f dx-\x({d-A I f + i's/7nl) R d-^, (1.39) 

- - A j 'dx-ti^ArAj} 2 ). (1.40) 



We can also construct the Noether charges corresponding to the supersymmetry trans- 



formation (|1.22| ). As in the three dimensional case it is convenient to decompose the 



supercharge in two components: 

Q + = PlQ, Q- = PrQ. 

The resulting eight component supercharges are given by 

Q + = 2 J dx'ti (#V fl 0_ A 7 ) , (1.41) 

Q- = -2gJ dx-ti (j + ^R + '-[AjAjWS - f3jf3f)^ R ) . (1.42) 

Finally we make a short comment on dimensional reduction of SYM 3+ i and SYM 5+ i. 
These systems can be constructed repeating the procedure just described. However there 
is an easier way to construct the Hamiltonian and supercharges for the dimensionally 
reduced theories, namely one has to truncate the unwanted degrees of freedom in the ten 
dimensional expressions. This is especially easy for the bosonic coordinates: one simply 
considers indices / and J running from one to two (for D = 4) or to four (for D — 6). 
The fermionic truncation can also be performed by requiring the spinor ipR to be 2- 
or 4-component. Then the only problem is the choice of 2 x 2 or 4 x 4 beta matrices 
satisfying 

{P I ,P J } = 2S IJ , (1.43) 

that can be done easily. 

2 Zero Modes and Light Cone Vacuum. 



The results of this section are based on the paper JlO | 



2.1 Gauge Fixing in DLCQ 

We consider the supersymmetric Yang-Mills theory in 1+1 dimensions [Q which is de- 



scribed by the action Ql.ll 



S = j d 2 x tr (—F^F"" + l -D^D»<P + i^D^ - 2ig0 lh ^\ . (2.1) 

A convenient representation of the gamma matrices is 7 = a 2 , 7 1 = ia 1 and 7 s = ct 3 
where a a are the Pauli matrices. In this representation the Majorana spinor is real. We 
use the matrix notation for SU(N) so that Afj and \IZjj are N x N traceless matrices. 

We now introduce the light-cone coordinates x ± = -j^(x° ± a; 1 )- The longitudinal 
coordinate x~ is compactified on a finite interval x~ G [-L, L] |63| , |68fl and we impose 
periodic boundary conditions on all fields to ensure unbroken supersymmetry. 

10 



The light-cone gauge A + = can not be used in a finite compactification radius, 
but the modified condition d^A + = |j^| is consistent with the light-like compactifica- 



tion. We can make a global rotation in color space so that the zero mode is diagonalized 
Afj(x + ) = Vi(x + )5ij with Sjfj = |52| . The gauge zero modes correspond to a (quan- 



tized) color electric flux loops around the compactified space. The modified light-cone 
gauge is not a complete gauge fixing. We still have large gauge transformations preserv- 



ing the gauge condition d-A + = 0. There are two types of such transformations [58], [59 
displacements Tjj and central conjugations Tq. Their actions on the physical fields of the 
theory and complete gauge fixing will be discussed in the end of this subsection. Now we 
just mention that being discrete transformations, To and Tc don't affect quantization 
procedure. 

The quantization in the light-cone gauge with or without dynamical A + is widely 
explored in the literature [71, 65, [M|, || [25| , here we provide only the results which are 



useful for later purposes. The quantization proceeds in two steps. First, we must resolve 
the constraints to eliminate the redundant degrees of freedom. There are two constraints 
in the theory, 

-D 2 _A~ = gj + , (2.2) 

y/2iD. X = g[Ml (2.3) 

where \1/ = (ip, x) T an d the current operator is 

J + (x) = i[0(x),D_0(x)] - -U^0r),V(x)}. (2.4) 

i V2 



Different components of (|2.2|) , Q2.3|) play different roles in the theory. First we look 
at diagonal zero modes of these equations. The diagonal zero mode of ( |2.3| ) gives us 
constraints on the physical fields: 

[<M°=0- (2.5) 

There is no sum over i in above expression. As one can see this constraint leads to 

o 
decoupling of Xu, this field plays the role of Lagrange multiplier for above condition. 

o ~ o 

The same is true for An, the corresponding constraint is Ju= 0. The reason we treated 



the diagonal zero modes of ( |2.2|) and (|2.3| ) separately is that for all other modes the 
D- operator is invertible and instead of constraints on physical fields ip and <p one gets 
expressions for non-dynamical ones: 

A ' = -W J+ ' * = s355-l**l- (2 ' 6) 

The next step is to derive the commutation relations for the physical degrees of 
freedom. As in the ordinary quantum mechanics, the zero mode Vi has a conjugate 
momentum p = 2Ld + Vi and the commutation relation is [t^Pj] = iSij. The off-diagonal 
components of the scalar field are complex valued operators with <pij = (</)ji)'. The 

11 



canonical momentum conjugate to 



'ij IS 71 {j 



(D_c 



iji- 



They satisfy the canonical 



commutation relations 71, p4| 



i i 

[(f)ij(x),TT kl (y)} x + =y + = [cf) ij (x),D^(j)i k (y)] x + =y + = -{5 ik 8ji - —5 ij 5 k i)5(x -y ). (2.7) 

On the other hand, the quantization of the diagonal component <fia needs care. As 



mentioned in [71], the zero mode of </>„, the mode independent of x~ 1 is not an independent 
degree of freedom but obeys a certain constrained equation pB| , 71, ^T| . Except the zero 
mode, the commutation relation is canonical 



{x),d.<f )jj {y)} x+=y+ = -(1 - —)8ij 



S(x 



r) ~^L 



(2.8) 



The commutator of diagonal and non-diagonal elements of <fi vanishes. The canonical 
anti-commutation relations for fermion fields are [^4 



1 1 



(2.9) 



There are two differences between this expression and one from |64| . First one is technical: 
we consider commutators for SU(N) group, this gives 1/N term. Second difference is 



that unlike |64| we include zero modes in the expansion of ip, we also include such modes 
in non-diagonal elements of 4>. 

Finally we return to the problem of complete gauge fixing. The actions of To and Tc 
on physical fields are given by |58|, [5!| : 



T, 



riiTC 



D 






E 



11; 



0. 



'c 



Vi(x~ 



,-ni(rii — rij)x 
exp { ^ ) ipij , 



L 



exp 



.myrii — rij)x 



L 



i.j- 



Vi'K 



Vi(x + ) H -, Vi = n(l/N - 5 iN ) 



(2.10) 



(2.11) 



tpij -> exp( 



m(vi 



Uj)X 



L 



■)V- 



u • 



exp 



■ KlXVi — Vj)X 



L 



There are also permutations of the color basis % —*■ P(i) which leave the theory invariant. 
These symmetries preserve the gauge condition d-A + = 0, but two configurations related 
by To, Tq or P are equivalent. To fix the gauge completely one therefore considers t>j 
only in the fundamental domain, other regions related with this domain by To, Tc or P 
give gauge "copies" of it |37|- The easiest thing to do is to describe the boundaries of 



fundamental domain imposed by displacements To: —j^ < v « < t^. The invariance 
under Tc limits this region even more, but since we will not need the explicit form of 
fundamental domain, we do not discuss such limits for SU(N). For the simplest case of 
SU{2) the fundamental domain is given by < v i - —v 2 < yj;, the result for SU(3) 



can be found in |59" j. The P symmetries do not respect the fundamental domain, so they 
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are not symmetries of gauge fixed theory. However there is one special transformation 
among P which being accompanied with combination of To and Tq leaves fundamental 
domain invariant. Namely if R is cyclic permutation of color indexes then there exists a 
combination T of Td and Tc such that S = TR is the symmetry of gauge fixed theory. 
The explicit form of T depends on the rank of the group, for SU(2) and SU(3) it may be 
found in |59[. The operator S satisfies the condition S N = 1 and it was used in classifying 
the vacua |59|, |70| . 



2.2 Current Operators 

The resolution of the Gauss-law constraint (|2.2| ) is a necessary step for obtaining the 
light-cone Hamiltonian. The expression for the current operator is, however, ill-defined 
unless an appropriate definition is specified, since the operator products are defined at the 
same point. We shall use the point-splitting regularization which respects the symmetry 
of the theory under the large gauge transformation. 

To simplify notation it is convenient to introduce the dimensionless variables z\ = 
Lgvi/ii instead of quantum mechanical coordinates v t describing A + . The mode-expanded 
fields at the light-cone time x + = are 



\j{X) 



J2 ay(7i)wy(n)e iknX + J2 a] i (n)u ij (-n)e ik " x ' , i^j, 




ik n x j_ y^ \ ( ri \ n ,.,(_ lrl \JknX- 
n=l 



<f>ii(x) = -=J2~r= [an(n)e tknX ~ + aUn)e ikn 
X 4tt ^ Vn v 

-1 / oo oo \ 

Mx) = - r7 =[J2bi j (n)e- ikn ^ + J2 b U n > iknX )^ ( 2 - 12 ) 



n=l / 

where k n = nir/L, Uij(n) = \/J\n — Zi + Zj\ []. The (anti) commutation relations for 
Fourier modes are found in [JTTJ, ^5| and in our notation they take the form 

[aij(n), 4/( m )] = sgn(n + Zj - Zi)S njm (Si k Sji - —S^m), 

{bij(n), b ] kl (m)} = 8 n>m {8 ik 8ji - —bijbki) (2.13) 

The zero modes in above relations deserve special consideration. Although we formally 
wrote them as eiy(O) and 6^(0), these modes also act as creation operators because the 
conjugation of zero mode gives another zero mode: 

4(0) = 0^(0), 6| 7 -(0) = 6 i< (0). (2.14) 

In particular the diagonal components of fermionic zero mode are real and we will use 
them later to describe the degeneracy of vacua. Now we concentrate our attention on 

Uij(n) is well-defined in the fundamental domain. Similarly, (-D-) 2 in the Gauss-law constraint have 
no zero modes in this domain. 
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non-diagonal zero modes. In the fundamental domain all z% are different, then one can 
always make take them to satisfy the inequality z^ < Zn-i < ■ ■ ■ < %\ in this domain. 
Such condition together with ( 2.13|) leads to interpretation of Oij(0) as creation operator 
if i < j and as annihilation operator otherwise. The situation for fermions is more 
ambiguous. One can consider 6y(0) as creation operator either when i < j or when 
% > j, both assumptions are consistent with Q2.13| ). Later we will explore each of these 
situations. 

Let us now discuss the definition of singular operator products in the current ( |2.4|) . 
We define the current operator by point splitting: 



J + = lim (j + ^(x; e) + J + ^(x; e)) 
where the divided pieces are given by 

J + *(x; t) -- 
J + ^(x;e) -- 



e- i ^ M 6(x 



(x 






V2 



(2.15) 

(2.16) 
(2.17) 



Here M is diagonal matrix: M = diag(zi, . . . , zn). An advantage of this regularization 
is that the current transforms covariantly under the large gauge transformation. 

To evaluate ( |2.16| ) and Q2.17| ) we will generalize the approach used in JZTJ, |65| to the 
SU(N) case. First let us calculate the vacuum average of bosonic current. Taking into 
account the interpretation of zero modes as creation-annihilation operators we obtain: 



(0|J± (x;e)|0> = : 



-e)D 

2L i z i- z k) . 



* IT \^"i ^J 



ik\X 

x~)ik\fy 



-e)D. 



x 



)kj 



AL 



-i — 

e 2^ 



k m>0 



'■'-•'■Me-*» e (5y-^- fc ) + 



k<j k>i 

Evaluating the sum and taking the limit one finds: 

1 



lim Jj- ,(x;e) 



-■■^ j ^)--+j L (^-(N + l-2t))5 ij , 



(2.18) 



(2.19) 



where : J + $ : is the naive normal ordered currents. To be more precise, we have omitted 
the zero modes of the diagonal color sectors in which the notorious constrained zero mode 
p3| appears. 

The result for fermionic current depends on our interpretation of zero modes as 
creation-annihilation operators and it is given by 



lim Jf; (x: e) 



-■■Jt J ^)---l L ^T(N + l-2i))5 ir 



(2.20) 
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The minus sign here corresponds to the case where %(0) is a creation operator if i < j 
(i.e. the convention is the same as for the bosons) and plus corresponds to the opposite 
situation. As can be seen, J + ^ and J + ^ acquire extra z dependent terms, so called 
gauge corrections. Integrating these charges over x~ , one finds that the charges are time 
dependent. Of course this is an unacceptable situation, and implies the need to impose 
special conditions to single out 'physical states' to form a sensible theory. The important 
simplification of the supersymmetric model is that these time dependent terms cancel, 
and the full current ( 2.15Q becomes 



W =: J $* : + : J t^ : +C ^r (2-21) 

Depending on the convention for fermionic zero modes the z independent constants Cj 
either vanish or they are given by 

Ci = ~(N + l-2i). (2.22) 

The regularized current is thus equivalent to the naive normal ordered current up to an 
irrelevant constant. Similarly, one can show that P + picks up gauge correction when the 
adjoint scalar or adjoint fermion are considered separately but in the supersymmetric 
theory it is nothing more than the expected normal ordered contribution of the matter 
fields. 

In one sense these results are a consequence of the well known fact that the normal 
ordering constants in a supersymmetric theory cancel between fermion and boson con- 
tributions. The important point here is that these normal ordered constants are not 
actually constants, but rather quantum mechanical degrees of freedom. It is therefore 
not obvious that they should cancel. Of course, this property profoundly effects the 
dynamics of the theory. 

2.3 Vacuum Energy 

The wave function of the vacuum state for the supersymmetric Yang-Mills theory in 1+1 
dimensions has already been discussed in the equal-time formulation [p 1 ?] . An effective 
potential is computed in a weak coupling region as a function of the gauge zero mode by 
using the adiabatic approximation. Here we analyze the vacuum structure of the same 
theory in the context of the DLCQ formulation. 

The presence of zero modes renders the light-cone vacuum quite nontrivial, but the 
advantage of the light-cone quantization becomes evident: the ground state is the Fock 
vacuum for a fixed gauge zero mode and therefore our ground state may be written in 
the tensor product form 

|fi) = $[z] ® |0>, (2.23) 

where we have taken the Schrodinger representation for the quantum mechanical degree 
of freedom z which is defined in the fundamental domain. In contrast, to find the ground 
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state of the fermion and boson for a fixed value of the gauge zero mode turns out to be 



a highly nontrivial task in the equal-time formulation ||67|| . 

Our next task is to derive an effective Hamiltonian acting on $[z\. The light-cone 
Hamiltonian H = P~ is obtained from energy momentum tensors, or through the canon- 
ical procedure: 

rr 9 2 L 1 ^ 8 . ,_d_ 

+ jy x -tr[- 9 ^J + ^J + + 0=[<p,^i<i ) ^] ] j, (2.24) 

K(z) = IIsm 2 (— ^4 (2-25) 

i>j 

where the first term is the kinetic energy of the gauge zero mode, and in the second 
term the zero modes of £)_ are understood to be removed. Note that the kinetic term 
of the gauge zero mode is not the standard form —d 2 /dz 2 but acquires a nontrivial 
Jacobian K which is nothing but the Haar measure of SU(N). The Jacobian originates 
from the unitary transformation of the variable from A + to v, and can be derived by 
explicit evaluation of a functional determinant |)B], |SD|. In the present context it is found 
in |nj]. Also we mention that Hamiltonian ( |2.24|) seems to contain terms quadratic in 



o 
diagonal zero modes ip u . However using constraint equations one can show that the total 

contribution of all such term vanishes. This also can be seen by using the fact that 

Hamiltonian is proportional to the square of supercharge (|2.34|) . 



Projecting the light-cone Hamiltonian onto the Fock vacuum sector we obtain the 
quantum mechanical Hamiltonian 

^-0^£|-*-w|- + K,, + ^, (2.26) 

where the reduced potentials are defined by 



Vjj = - y - J dx~{trJ + - T J + ), (2.27) 



g 2 f L i 

V* ^^f^dx-^M^M), (2-28) 

respectively. As stated in the previous subsection, the gauge invariantly regularized 
current turns out to be precisely the normal ordered current in the absence of the zero 
modes. It is now straightforward to evaluate Vjj and V^p in terms of modes. One finds 
that they cancel among themselves as expected from the supersymmetry: 



Vjj — —V<t>ip 



16tt 2 



y y y — + 

n^Ll T^( n - Z i + Z k) ( m + Z 3 ~ Z k) n ^Lx «™ 
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°° I 1 

hi y vfc ^ ~ 2i + Zk ^ Zj ~ Zk) hi ( n + z i - z k) ( z k -Zi)' 

ij i>k>j \ Zk ~ Z i)\ Z J ~ Z k) 



(2.29) 



This cancellation was found as the result of formal manipulations with divergent series 
like ones in the right hand side of the last formula. Such transformations are not well 
defined mathematically and as the result they may lead to the finite "anomalous" con- 
tribution. The famous chiral anomaly initially was found as the result of careful analysis 
of transformations analogous to ones we just performed (T|. However if one considers 
derivatives of Vjj or V^ with respect to any Zi then all the sums become convergent, 
the order of summations becomes interchangeable and as the result the derivatives of 
Vjj + V^ vanish. Thus if there is any anomaly in the expression above it is given by z- 
independent constant. Such constant in the Hamiltonian would correspond to the shift 
of energy levels and usually it is ignored. However in supersymmetric case there is a 
natural choice for such constant: in order for vacuum to be supersymmetric it should be 
zero. Below we assume that SUSY is not broken, then we expect that ( |2.29| ) is true. 
Thus we arrive at 

The relevant solutions of this equation should be finite in the fundamental domain, this 
requirement leads to discrete spectrum due to the fact that Jacobian vanishes on the 
boundary of this domain. However the operator H Q is elliptic, and therefore it can't have 
negative eigenvalues. If the eigenvalue problem 

H $(z) = E®{z) (2.31) 

has a solution for E — 0, this solution corresponds to the ground state of the theory. 
It is easy to see that such solution exists and it is given by $(z) = const 0. We have 
thus found that the ground state has a vanishing vacuum energy, suggesting that the 
supersymmetry is not broken spontaneously. 

2.4 Supersymmetry and Degenerate Vacua. 

As we saw in the previous subsection supersymmetry leads to the cancellation of the 
anomaly terms in current operator. However these terms played an important role in 
the description of Z^ degeneracy of vacua 175] , so we should find another explanation of 



this fact here. It appears that fermionic zero modes give a natural framework for such 
treatment. 



2 some authors prefer to rewrite this to include the measure in the definition of the wave function and 
then in SU(2) for example the ground state wave function is a sin 
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First we will generalize the supersymmetry transformation given in |64j to the present 
case, i.e. we include A + and the zero modes of fermions. The naive SUSY transformations 
spoil the gauge fixing condition, so we combine them with compensating gauge transfor- 
mation following [Q. In three dimensional notation (spinors have two components and 
indices go from to 2) the result reads: 

5 A, = Igfy* - D„Ig7_ J-*, (2.32) 



The difference between above expression and those in [64] is that we include the zero 
modes. Namely we defined \I> as the complete field with all the zero modes included and 
\l/ as fermion without diagonal zero modes. The introducing of \& is necessary, because 
diagonal zero modes form the kernel of operator D_, so -^- is not defined on this subspace. 
In particular we are interested in supersymmetry transformations for A + and fermionic 
zero modes. Performing a mode expansion one can check that diagonal elements of matrix 
[jj-ip,ip}° vanish, then from (|2.32j) we get: 



Ht = 


% T o 

7f + ^ 





-2d + A+e + . 



(2.33) 

This expression is written in two component notation and the decomposition of spinor 
e: e = (e+, S-) T is used. Note that since eQ = V2(e + Q~ + £-Q + ) the fields involved in 
transformations ( p.33| ) don't contribute to Q + , this is consistent to the fact that being 
x~ independent they don't contribute to P + . The equations ( |2.33| ) look like supersym- 



metry transformation for the quantum mechanical system built from free bosons and free 
fermions. In fact as one can see the supercharge Q~ is the sum of supercharge for the 
quantum mechanical system and from the QFT without diagonal zero modes: 

Q~ = -2g f dx~tr(J + —4>) + ALtr(d + A + ip). (2.34) 

Calculating {Q~) 2 and writing the momentum conjugate to A + as differential operator 
we reproduce Hamiltonian ( |2.24| ). Note that ip there has all the zero modes in it. The 
square of another supercharge 

Q + = 2 I ' dx-ti(i)D_<p) (2.35) 

gives P + while the anti-commutator of Q~ with Q + is proportional to the constraint 
(B3I) and thus vanishes. 



3 using Schrodinger coordinate representation for quantum mechanical degree of freedom - note that 
the QFT term has non-trivial dependence on the quantum mechanical coordinate. 

18 





One can check that although [-0^, H] does not vanish, this commutator annihilates 

Fock vacuum |0), then it also annihilates \Q). In subsection 1 we mentioned that Xu 

decouples from the theory, and therefore it commutes with Hamiltonian. Thus acting on 

o o 
the vacuum state \fl) by diagonal elements of either ip or X we get states annihilated by 

P~ and P + (the latter statement is obvious since zero modes commute with momentum). 

Not all such states however may be considered as vacua. Although we fixed the gauge in 

subsection 1, the theory still has residual symmetry P, corresponding to permutations 

of the color basis. Physical states are constructed from operator acting on the physical 

vacuum \Q) and both the operators and the physical vacuum must be invariant under 

P. Such objects can always be written as combinations of traces. The candidates for the 

o o 

vacuum state may have any combination of ip and X inside the trace, here and below we 

o 
consider only diagonal components of zero modes. Since X is not dynamical we have the 

usual c-number relation 

{Xu,X 3J } = (2.36) 





instead of canonical anti-commutator, so XX= 0. From the relations ([2.1 3|) one finds: 



^izW 1 -^' (2 ' 37) 

o o o o o 

also we have Xtp= — ipX. Using all these relations and the SU(N) conditions tr(0>) = 

o o o 

and tr(x) = we find that the only nontrivial trace involving only zero modes is tr(0>X). 

Then the family of vacua is given by: 

( °o \ n 

tr(^x) \D), 0<n<N-l. (2.38) 

o 
The region for n is determined taking into account the fact that X is anti-commuting field 

with N — 1 independent components. Thus we explained the Z N degeneracy of vacua 

first mentioned in ||81|| . 

In addition to this discrete vacuum degeneracy supersymmetric theories also have a 

continuum space of vacua called moduli space. In DLCQ approach the moduli space 

is easy to understand. Let us suppose that scalar field developed a VEV. To have 

a consistent theory this VEV should commute with the Wilson loop in the compact 

direction, which in our case happened to be exp(i f dx~A + ). Since A + is a general 

diagonal matrix this leads to the condition for the VEV: (<frij) = WiSij. Now we can make 

the substitution — > + (0) in the supercharges ( p.34|) to find the correction in Q~ due 

to the scalar VEV: 

5Q- = -2igti (w f dx~[(j)ip]\ . (2.39) 

We used integration by part and the equation D_w = 0. Taking into account the 
fermionic constraint ( |2.5| ) we conclude that for any diagonal w. 5Q~ = 0, i.e. we can 
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choose the state with arbitrary VEV (<pij) as the new vacuum. This is precisely the 
moduli space of the theory: the models constructed starting from different vacua are not 
coupled with each other. 

2.5 Solving for Massive Bound States. 

As we saw the zero modes play an important role in the description of the vacuum. 
However solving for massive bound states one usually neglect the zero mode contribution. 
Does this lead to errors in the mass spectrum? The answer depends on the problem we are 
solving. If one is interested in the spectrum of the theory at the finite value of resolution 
then zero modes are important, but as we will show their contribution becomes smaller 
and smaller as the resolution goes to infinity, so they may be neglected if one is interested 
only in the large K extrapolation. 

First let us formulate the DLCQ problem with zero modes precisely. We will use the 
Hamiltonian formulation, but the consideration for SDLCQ formalism is the same. The 
space of states of the theory is the direct product of Fock space and quantum mechanical 
Hilbert space for zero modes: a general state can be written as 

\state) = $(z) <g> \FockState), (2.40) 

the Hamiltonian has the form: 

H = K(z)+V(z,a,a\b,tf). (2.41) 

Here K(z) is some differential operator, while V is some function of zero modes z and 
creation-annihilation operators (see for example Q2.24J) ). In general one should solve the 



bound state problem H\ty) = E\^f) in two steps: first one should determine the effective 
potential V: 

V(z,a,a\b,b j )\FockState) = V(z)\\Fock State) (2.42) 

and then solve the Schrodinger equation for zero modes: 

(K(z)+V(z))V(z) = EV(z). (2.43) 

However in practice this is hard to carry out. Fortunately, solving the Schrodinger 
equation is not important for calculating the continuum limit of mass spectrum. The 
reason for this is the following. 

Studying the large L limit in DLCQ one is usually interested in situation when the 
total momentum P + = J2 n i/L is kept fixed. Then most of the terms in V (and thus in 
V) are of order L°, while K{z) scales like L. Assume for a moment that the whole V is 
of order one, then one can consider V as perturbation and use the standard expression 
for the eigenvalue: 

E % = Ef ] + J dz*l(z)V(z)*i(z), ( 2 - 44 ) 
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where E\ and ^i(z) are eigenvalue and eigenfunction of unperturbed system. To get 
finite masses in the continuum limit only the ground state of K{z) should be considered: 
i = and Eq = in the last expression. Introducing the averaging procedure as 

(A) = J dz*t(z)A(z)Mz) ( 2 - 45 ) 

we find: E = (V) and thus the continuum eigenvalues are just solutions of the z- 
independent equation: 

(V{a, a\ b, b ] ))\FockState) = E\FockState). (2.46) 

The assumption of L° scaling for V is not the trivial one. Namely it is responsible for 
the difference in the constraint equations in DLCQ and continuum cases. For example 



looking at the Hamiltonian Q2.24J ) one can see that V(z) includes a term linear in L: 



^^^4(0)^(0), (2.47) 

so the assumption being false for V may be satisfied for V only dynamically. One 
can make this specific term vanish if instead of DLCQ constraint f dx~Ju(x) = its 
continuum version 

fdx-Jij(x) = (2.48) 

is used. Of course imposing this condition is not enough to make all the terms in V 
to be of order L°, but following the usual path in DLCQ calculations we choose not to 
impose other conditions explicitly. In our numerical study we rather perform calculations 
with Hamiltonian (V(a, a\ b, b')) in the sector satisfying (|2.48| ) and then concentrate our 
attention only on states whose masses can be extrapolated to finite value. This way we 
make sure that our assumption V ~ 1 holds and thus the z dependence is not important. 
To summarize, we have shown that zero modes of gauge field and diagonal zero modes 
of fermions play an important role in the description of vacuum structure. However if 
studying the bound state problem for the states with nonzero total momentum P + one 
is interested only in the extrapolation to the continuum limit, the zero mode of A + can 
be omitted from the theory. This fact leads to significant simplifications in the numerical 
procedure. As soon as A + is excluded from the theory one also has to exclude the bosonic 
zero modes (otherwise the expression 1/0 is encountered in the ( 2.12Q ). What about the 



fermionic zero modes? In principle we can either keep them or disregard them. However 
in the latter case one should be very careful: as we will see in the next section such modes 
play an important role in the ensuring of supersymmetry. 

3 Fermionic Zero Modes and Exact Supersymmetry. 

In this section we study the relation between conventional DLCQ and its supersymmetric 
version. Since usual DLCQ is formulated for the Hamiltonian we should rewrite SDLCQ 
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in the same form. Here one encounters the first difference between two schemes: in 
DLCQ the fermions can be chosen to be either periodic or antiperiodic on x~ , but in the 
Hamiltonian formulation of SDLCQ they must be periodic due to supersymmetry. Then 
one encounters the problem of fermionic zero modes. However the boundary conditions 
is not the only difference between the two approaches. Even after we choose periodic 
fermions, DLCQ still has an ambiguity emerging from the choice of regularization scheme. 
Taking the simplest SUSY system as an example we will show that supersymmetry 
dictates the unique regularization and we study the relation between this prescription 
and the principal value scheme, which is usually used in the DLCQ calculations. We 
show that fermionic zero modes play an important role in deriving this relation. 

As we already mentioned in section one the simplest supersymmetric system in two 
dimension is the one involving only gauge fields and adjoint fermions ||57| . We derive all 
the relations for this particular system. 

3.1 Zero Modes and Supersymmetric Regularization. 

We consider the 1 + 1 dimensional SU(iV) gauge theory coupled to an adjoint Majorana 
fermion. The light-cone quantization of this model in the light-cone gauge and large iV 



limit has been dealt with explicitly before [^9[ 20f| . The expressions for the light-cone 



momentum P + and light-cone Hamiltonian P for this model are 

P + = [ dx-tr(iV2^ij), (3.1) 



Here J^ = —\<2tpikipkj is the longitudinal current. It is well known that at a special 
value of fermionic mass (namely m 2 = g 2 N/n) this system is supersymmetric [^j. This 
special value of the fermion mass will be denoted by rnsusY- At this supersymmetric 
point, the supercharge is given by 

Q- = V2gj dx-tr(ijij—^) (3.3) 

which satisfies the supersymmetry relation {Q~,Q~} = 2\/2P~. This may be checked 
explicitly by using the anticommutator at equal x + : 

{i>ij{x~),^ki{y~)} = -S u 8 jk 8(x- - y~). (3.4) 

In the DLCQ formulation, the theory is regularized by a light-like compactification, and 
either periodic or antiperiodic boundary conditions may be imposed for fermions. If 
P + denotes the total light-cone momentum, light-like compactification is equivalent to 
restricting the light-cone momentum of partons to be non-negative integer multiples of 
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P + /K, where K is some positive integer that is sent to infinity in the decompactified 
limitQ. Anti-periodic boundary conditions will in general explicitly break the supersym- 
metry in the discretized theory, although supersymmetry will be restored in the decom- 
pactification limit K — > oo |20f . If we wish to maintain supersymmetry at any finite K, 



we must at least impose periodic boundary conditions for the fermions. This, however, 
leads to the notorious "zero-mode problem"^ From a numerical perspective, omitting 
zero-momentum modes in our analysis is absolutely necessary, since it guarantees a finite 
Fock basis for each finite resolution K. The mass spectrum of the continuum theory 
may be obtain by extrapolating from a sequence of finite mass matrices M 2 = 2P + P~ . 
But are we really justified in omitting the zero- momentum modes? To date, the gen- 
eral consensus is that omitting zero momentum modes in a two dimensional interacting 
field theory does not affect the spectrum of the decompactified theory, where K — > oo. 
Actually, the numerical results of the next subsection are consistent with this viewpoint. 
However, the goal of this work is to understand the structure of a supersymmetric 
theory at finite resolution. As we will see shortly, understanding why the DLCQ and 
SDLCQ prescriptions differ involves studying certain intermediate zero-momentum pro- 
cesses. But first, we need to be more precise about the form of the light-cone operators 
of the theory. If we expand the fermion field ipij in terms of its Fourier components, we 
may express the uncompactified light-cone supercharge and Hamiltonian in a momentum 



space representation involving fermion creation and annihilation operators: (j57|, |29|, |20|): 

i2 -1 / 4 o r°° ( 1 1 1 \ 

Q~ = -=- / dk 1 dk 2 dk 3 5(k 1 + k 2 - k 3 ) — + -. — x 

y/n Jo \kx k 2 k 3 J 

x (4(^i)4#2)M^) + &l# 3 )M^i)M^)) , 

(3.5) 
m 2 f°° dk , g 2 N f°° dk f k k 



m? f°° dk lf .... ... g*l\ f°° dk f K , fc , t/n , ... 

g 2 f°° r t t 

+— - / dk 1 dk 2 dk 3 dk 4: 5{ki+k 2 -k 3 -k A )A(k)b\Ak 3 )b] i (k i )b k i(ki)b H (k 2 ) + 
Zn Jo L 

+5(k 1 +k 2 +k 3 -k 4 )B(k) x 

x (bUkJhiih^ui^ih) -bljih^ik^bU^buih)] 



(3.6) 



with 



A(k) 



(h - k 2 ) 2 (h + k 2 ) 



B ^ = WTW-WTW- (3 ' 7) 



4 K is sometimes called the 'harmonic resolution', or just 'resolution'. 

5 For anti-periodic boundary conditions, the light-cone momentum of partons is restricted to odd 
integer multiples of P + / ' K, and so there are no zero-momentum modes in such a formulation. 
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As we mentioned earlier, the continuum theory is supersymmetric for a special value 
of fermion mass. We will therefore consider only the case m = rnsusY- in the DLCQ 
formulation, one simply restricts integration of the light-cone momenta k{ in expression 
Q3.6| ) for P~ above to be positive integer multiples of P + /K. i.e. one simply drops the 
zero-momentum mode. The DLCQ mass spectrum is then obtained by diagonalizing 
the mass operator M 2 = 2P + P~. Similarly, in the SDLCQ formulation, the light-cone 



momenta k{ in expression (|3.5|) for Q~ are restricted to positive integer multiples of P + /K. 
One then simply defines P to be the square of the supercharge: 2y / 2-P~ = {Q~,Q~}- 
The mass operator M 2 = 2P + P~ is then easily constructed and diagonalized to obtain 
the SDLCQ spectrum. 

In general, the following observations are made; at finite resolution, the DLCQ spec- 
trum of a supersymmetric theory is not supersymmetric. However, supersymmetry is 



restored after extrapolating to the continuum limit K — > oo (see [20], for example). In 
contrast, for any finite resolution, the SDLCQ spectrum is supersymmetric. The DLCQ 
and SDLCQ spectra agree only in the decompactified limit K —>■ oo. 

Not surprisingly, the difference in the DLCQ and SDLCQ prescriptions at finite reso- 
lution may be understood as a zero-mode contribution. What is surprising is that we can 
encode the effect of these zero-mode contributions into a simple well defined operator. 
The main result here is the the precise operator form of this contribution at finite K. 

In order to motivate our argument, note that the anticommutator for the supercharge 
Q~ in the continuum theory involves products of terms of the form h< (kp (0)b(k) and 
tf (p)b(0)b(p) , and these provide contributions to P~ that may be expressed in terms of 
non-zero momentum modes. The problem is exacerbated by the fact that the coefficients 
of these terms behave singularly. To examine this more closely, we consider the discretized 
theory where the light-cone momenta hi in the expression for Q~ [eqn( p.5| )] are restricted 
to positive integer multiples of P + /K. We also include the effects of zero-momentum 
modes by introducing an 'e regulated zero mode', which are modes with momentum 
hi = e, where e is much less than P + /K, and is sent to zero at the end of the calculation. 
Then the anticommutator of two Q~ gives contributions of the following form: 

)b\k)b\e)b(k + e), (- + e p(p + e)b(e)b(p) 1 = 




k(k + e) e p(p + e 

tf{k)b(k+ep(p+e)b(p) 



11 1 e 2 

+ 



c 2 



p(p+e) k(k+e) pk(p+e)(k+e) 



(3.8) 



where any terms involving an e regularized zero mode on the right-hand-side are dropped 
and zero modes are omitted from P~ . We have suppressed all matrix indices in this 
expression. In the limit e — ^ the last term on the right-hand-side in the brackets 
vanishes, while the first term is the pure momentum-independent divergence that was 
identified in an earlier study of this model p0 |, and is canceled if we adopt a principal 



value prescription for singular amplitudes in the definition of P . The second term 
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however, is clearly a finite contribution to P~, although it arises from the e regulated 
zero modes in Q~ , which are not present in the SDLCQ prescription for defining Q~ . 
Consequently, in order to ensure the supersymmetry relation {Q~ , Q~} = 2\J2P~ in the 
discretized formulation, we must include an e regularization of the zero modes in the 
definition for Q~ , and then apply a principal value prescription in the presence of any 
singular processes to eliminate 1/e divergences. 

Stated slightly differently, we may decompose the supercharge into a part without 
zero modes Qsdlcq 0- e - ^* = nP + /K,n = 1,2, . . .), and terms with e regularized zero 
modes, Q~ . The anti-commutator {Qsdlcq^Q^} contains only terms with e regulated 
zero-modes. Since Q~ = Qsdlcq + Qe one finds 

\Q SDLCQ' Q SDLCQ J = ^v2P SDLC q = ^V^PdlCQ ~~ \Q e 'Qe \ PV i (3-9) 

after dropping any e regulated zero- mode terms in the calculated expression for {Q~ , Q~}. 
Note that the first equality above is just the definition for the light-cone Hamiltonian 
P in the SDLCQ prescription. The PV abbreviation on the right hand side indicates 
a principal value regularization prescription, which is tantamount to dropping all 1/e 
terms as e — > 0. The procedure is well known in the context of the present model 
|p0| . It is clear that our definition for Psdlcq gives rise to the supersymmetry relation 
[Q~sdlcqi Psdlcq} = 0' which yields a supersymmetric spectrum for any finite resolution 
K. Moreover, we know that Psdlcq anc ^ ^dlcq yield the same spectrum in the contin- 
uum limit K — > oo, so it remains to calculate the difference at finite resolution K. We 
will write this difference in terms of their respective mass operators: M 2 = 2P + P~ . A 
straightforward calculation of the anticommutator on the right-hand-side of ( |3.9| ) leads 
to the result: 

MIdlcq ~ M 2 DLCQ = Ml = - 9 -^- £ i4(n)B<,-(n) 

_g_ 2 NK J_ + ^ B l {m)B t. in)Bklim)Bli{n y (3.10) 

We also write down the expression for M^ LC q in the theory with periodic fermions: 

Mdlcq = ^ E BUn)B l3 (n)^ + £ —^) + 
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(n 2 -n A ) 2 (n 1 +n 2 ] 
1 1 



it r„_\rt 
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Sl i (n3)S} 4 (n4)S«(n 1 )S«(n 2 ) 

x (3.11) 



(^■MBmM^M^M - Bl i in 1 )B] l {n 2 )Bl{n & )B hi 

2 

In this expression the variable x = ZI ^ is a dimensionless mass parameter, and for the 
supersymmetric point we have x = 1. The sums are performed over positive integers, 
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< rii < K, and we employ a principal value prescription in sums labeled as YJ', which 
implies that terms of the form l/(k — k) 2 are dropped. In the SDLCQ procedure we 
calculate Q~ which is non-singular and requires no principal value prescription. 

The term M\ appears to be non-trivial due to the presence of B^B^BB terms on 
the right hand side of ( |3.10| ). However, the action of this term on any SU(iV) Fock state 
turns out to be equivalent to the first term, although with opposite sign, and twice the 
magnitude. Thus the action of the right hand side of (|3.10|) is equivalent to the single 
quadratic operator: 

Ml = ^ £ U(n)B(,(n). (3.12) 



Fortunately, we are able to test this analytical result by performing direct numerical sim- 
ulations of this model using both prescriptions, and comparing the differences observed 
with the above prediction. Interestingly, although this result was derived for large N, 
agreement turns out to be perfect for both finite and large N, which was verified using 



the finite N DLCQ algorithms developed in ||13| . We discuss this further in the next 
subsection. 

3.2 Soft SUSY Breaking and Numerical Results. 

In this subsection we compare the numerical results for different regularization schemes. 
Although in the continuum limit both PV and SUSY prescription ns should give the same 
results, the convergence of the masses as K — > oo might be different. So if at a given value 
of K one wants to get a better approximation to continuum masses, one scheme might 
work better than other. In previous subsection we described two regularization schemes 
and found the operator M\ describing the difference between them. It is convenient to 
introduce the family of regularizations labeled by parameter Y: 

P y =Ppv + YMI (3.13) 

Then at two special values of Y we get the PV and SUSY prescriptions: Pp V = -Py =0 , 
Psusy = Py=\- Since Pp V is defined for arbitrary value of fermionic mass m (not only for 
super symmetric one) the last equation also defines the family of regularizations beyond 
the SUSY point. On the other hand shifting the fermionic mass from its supersymmetric 
value is equivalent to introducing additional fermionic mass, i.e. to the soft SUSY break- 
ing. Below we will give a numerical results for bound state masses in the theory with 

2 

two new " coupling constants" : X = I ^ and Y which determine the value of fermionic 
mass and regularization scheme accordingly. 

Our investigation of this theory indicates that at X = 1 ( the supersymmetric value 
of the fermion mass) the lightest fermionic and bosonic bound states are degenerate with 
continuum masses approximately M 2 = 26 |2T| |J. Using P^usy we arrive at the same 
conclusion for any value of Y . 



Boorstein and Kutasov 2^ have investigated 'soft' supersymmetry breaking for small 



values of this difference, X — 1 and they found that the degeneracy between the fermion 
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and boson bound state masses is broken according to 

M 2 F (X) - M 2 B (X) = (1 - X)M B {1) + 0{(X - l) 3 ). 



(3.14) 



They calculated these masses using the PV prescription (Y = 0) with anti-periodic 
BC and found very good agreement with the theoretical prediction. We have compared 
this theoretical prediction at Y = 1 and we find that eq ( |3.14| ) is very well satisfied. At 
resolution K = 5, for example, the slope is 4.76 and the predicted slope Mb(1) is 4.76. 
The indication is that this result is true for any value of Y. 
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Figure 1: (a) The contour plots of Y — Y(X) for the mass squared of the lowest bound 
state in units of g 2 N/iT as a function of X = mir/g 2 N and Y (b)The contour plots of 
Y = Y(X) for the mass squared of the second lowest bound state in units of g 2 N/ir as a 
function of X = mm / ' g 2 N and Y (b) 

In Fig. 1 we show the contour plots of the mass squared M 2 of the two lightest 
bosonic bound states as a function of X and Y at resolution K = 10. These contours 
are lines of constant mass squared. Selecting a particular value of the mass of the first 
bound state then fixes a particular contour in Fig. la as a contour of fixed mass, which 
we can write as Y = Y p (X). 

Interestingly, constructing the same contour plot for the next to lightest bosonic 
bound state - see Fig. lb - yields contours that have approximately the same functional 
dependence implied by Fig. la. In fact, one obtains approximately the same contour plots 
for the next twenty bound states (which is as far as we checked). The simple conclusion 
is that the coupling Y which represents the strength of the additional operator affects all 
bound state masses more or less equally. This in turn suggests that at finite resolution, 
we can smoothly interpolate between different values of fermion mass X, and different 
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prescriptions specified by the coupling Y, without affecting too much the actual numerical 
spectrum. Of course, in the decompactification limit K — > oo, such a dependence on Y 
disappears, due to scheme independence. 

Since the lightest bosonic bound state is primarily a two particle state it is reasonable 
to truncate the Fock basis to two particle states. This will permit very high resolutions, 
which will be needed to carefully scrutinize any possible discrepancies between the two 
versions of 'soft' symmetry breaking presented here. In fact, we are able to study the 
theory for K up to 800. The mass of the lowest state as a function of the resolution for 
various values of X and Y are shown in Fig. 2. Each converging pair of lines - which 
extrapolate the actual data points - in Fig. 2 corresponds to different values of fermion 
mass X. The top upper curve in each pair runs through data points that were calculated 
via SDLCQ (i.e. Y — 1), while the lower corresponds to the PV (i.e. Y — 0) prescription 
commonly adopted in the literature. We find that each pair of curves converge to the 
same point at infinite resolution, although this may not be completely obvious for the 
lowest pair in the figure (corresponding to the critical mass X = 0). 

Away from X = 0, the SDLCQ formulation is fitted with a linear function of 1/K, 
while the PV formulation is fit with a polynomial of 1/K 2 ", where /3 is the solution of 
1 — X/2 = Trf3Cot(ir(3) jn|. It now appears that SDLCQ not only provides more rapid 
convergence for supersymmetric models, but also for the massive t'Hooft model, which 
is not supersymmetric. For the massless case, the situation is reversed; the SDLCQ 
formulation converges slower. It is fit by a polynomial in l/Log(K) and gives the same 
mass at infinite resolution as the PV formulation. This behavior may be understood 
from the observation that the wave function of this state does not vanish at x — 0. We 
have looked closely at 'small' masses, such as X = .1, and one finds that both PV and 
SDLCQ vary as a polynomial in 1/K 213 at large resolution. Thus careful extrapolation 
schemes must be adopted at small masses. 

We therefore conclude that the continuum of regularization schemes that interpo- 
late smoothly between the SDLCQ and PV prescriptions - which we characterized by 
the parameter Y - yield the same continuum bound state masses, although the rate of 
convergence of the DLCQ spectrum may be altered significantly. This implies that the 
contour plots observed in Fig. 1 eventually approach lines parallel to the Y axis, and the 
sole dependence on the parameter X is recovered. 

Interestingly, since the two-body equation studied here for the adjoint fermion model 
is simply the t'Hooft equation with a rescaling of coupling constant, we have arrived at 
an alternative prescription for regulating the Coulomb singularity in the massive t'Hooft 
model that improves the rate of convergence towards the actual continuum mass. Thus, 
a prescription that arises naturally in the study of supersymmetric theories is also appli- 
cable in the study of a theory without supersymmetry. We believe that this idea deserves 
to be exploited further in a wider context of theories. In particular, it is an open question 
whether this procedure could provide a sensible approach to regularizing softly broken 
gauge theories with bosonic degrees of freedom, and in higher dimensions. 

In any case, it appears that the special cancellations afforded by supersymmetry - 
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Figure 2: Mass of the of the lowest bound state in units of g 2 N/7r calculated in the 
t'Hooft model. The top pair is at X = 1, the second is at X = .5, and the bottom pair 
is at X = 

especially in the context of DLCQ bound state calculations - might have scope beyond 
the domain of supersymmetric field theory. This would be a crucial first step towards a 
serious non-perturbative study of theories with broken supersymmetry. 

4 Massless States in Two Dimensional Models. 

In this section we will study the structure of bound states for two dimensional super- 
symmetric models defined in section 1. We will concentrate most of the attention on 
the model obtained by dimensional reduction from SYM 2+ i. For this theory we will 
prove that any normalizable bound state in the continuum must include a contribution 
with arbitrarily large number of partons. By generalizing this proof to the theories with 
extended SUSY we show that this is the general property of supersymmetric matrix 
models. This scenario is to be contrasted with the simple bound states discovered in 
a number of 1 + 1 dimensional theories with complex fermions, such as the Schwinger 
model, the t'Hooft model, and a dimensionally reduced theory with complex adjoint 
fermions [T2|, |55[ . We also study the massless states of SYM 2+ i in DLCQ. Some of them 
are constructed explicitly and the general formula for the number of massless states as 
function of harmonic resolution is derived for the large N case. This section is based in 
part on the results of [§]. 
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4.1 Formulation of the bound state problem. 

The light-cone formulation of the supersymmetric matrix model obtained by dimension- 
ally reducing J\f = 1 SYM 2 +i to 1 + 1 dimensions was initially given in fMf , and it was 
summarized in the section 2 of these lectures. We simply note here that the light-cone 
Hamiltonian P~ is given in terms of the supercharge Q~ via the supersymmetry relation 
{Q-,Q-} = 2V2P-, where 

Q~ = 2 3 / 4 s J dx-ti |(i[0,a_0] + 2^)^-vl • (4- 1 ) 

In the above, 0y = 0jj(x + , x - ) and ip^ = tpij(x + ,x~) are N x N Hermitian matrix 
fields representing the physical boson and fermion degrees of freedom (respectively) of 
the theory, and are remnants of the physical transverse degrees of freedom of the original 
2 + 1 dimensional theory. This is a special feature of light-cone quantization in light- 
cone gauge: all unphysical degrees of freedom present in the original Lagrangian may be 
explicitly eliminated. There are no ghosts. 

In order to quantize and ip on the light-cone, we first introduce the following ex- 
pansions at fixed light-cone time x + = (the continuum counterpart of fl2.12|) : 



Mx-,0) = -±= r -tKf {^(k + )e-* k+x - + aUk + y k+x -) ; (4.2) 

V27T Jo \?2k + v ' 

^(x-, 0) = -L H dk + (b l3 (k + )e^ k+x - + b} i (k + )e ik+ ^) . (4.3) 

2\/TC Jo v J / 



2v% Jo 
We then specify the commutation relations 

>«(P + ), 4(9 + )] = {MP + )> 4(<Z + )} = <Hp + " ? + )M* (4-4) 

for the gauge group U(iV), or 

,(p + ), «L(<? + )] = {Mp + )> 4(<? + )} = <Kp + - <?+) (<W - ^M«) ( 4 - 5 ) 



Ojj \ 



for the gauge group SU(iV)f]. 

For the bound state eigen-problem 2P + P~\ty >= M 2 |\l/ >, we may restrict to the 
subspace of states with fixed light-cone momentum P + , on which P + is diagonal, and so 
the bound state problem is reduced to the diagonalization of the light-cone Hamiltonian 
P~ . Since P~ is proportional to the square of the supercharge Q~, any eigenstate |\& > 
of P~ with mass squared M 2 gives rise to a natural four-fold degeneracy in the spectrum 
because of the supersymmetry algebra — all four states below have the same mass: 

|tt>, Q + \^>, Q~\*>, Q + Q-\^>. (4.6) 



6 We assume the normalization tr[T a T b ] = S ab , where the T a 's are the generators of the Lie algebra 
ofSU(iV). 
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Although this four-fold degeneracy is realized in the continuum formulation of the the- 
ory, this property will not necessarily survive if we choose to discretize the theory in 
an arbitrary manner. However, a nice feature of SDLCQ is that it does preserve the 
supersymmetry (and hence the exact four-fold degeneracy) for any resolution. 

Focusing attention on zero mass eigenstates, we simply note that a massless eigen- 
state of P~ must also be annihilated by the supercharge Q~ , since P~ is proportional to 
(Q~) 2 - Thus the relevant eigen-equation is Q~\^ >= 0. We wish to study this equation. 
However, first we need to state the explicit equation for Q~ , in the momentum represen- 
tation, which is obtained by substituting the quantized field expressions (|4.2| ) and 
directly into the the definition of the supercharge ( |4.1| ). The result is: 

i2" 1 /4< 



Q~ = — -=— / dkidk 2 dk^8(ki + k 2 — k 3 
V 71 " J Q 
1 k 2 — ki 



4fc(fci)4i (fo)bij(fa) ~ b\j{h)aik{ki)a kj {k 2 )} 

(^ + y 2 - ^)P&(*i)H,-(k)M*s) + &1# 3 )M^)M^)] 



2 y /k 1 k 2 h 
1 h + k 3 

1 k 2 + k 3 



(4.7) 



In order to implement the DLCQ formulation [B3, K3| of the theory, we simply restrict 



the momenta fci, k 2 and k 3 appearing in the above equation to the following set of allowed 
momenta: {^-, x~, ^R~i . . .}. Here, K is some arbitrary positive integer, and must be 
sent to infinity if we wish to recover the continuum formulation of the theory. The 
integer K is called the harmonic resolution, and 1/K measures the coarseness of our 
discretization. Physically, 1/K represents the smallest unit of longitudinal momentum 
fraction allowed for each parton. As soon as we implement the DLCQ procedure, which 
is specified unambiguously by the harmonic resolution K, the integrals appearing in 
the definition of Q~ are replaced by finite sums, and the eigen-equation Q~\^ >= is 
reduced to a finite matrix problem. For sufficiently small values of K (in this case for 
K < 4) this eigen-problem may be solved analytically. For values K > 5, we may still 
compute the DLCQ supercharge analytically as a function of N, but the diagonalization 
procedure must be performed numerically. 

For now, we concentrate on the structure of the zero mass eigenstates for the contin- 
uum theory. Firstly, note that for the U(iV) bound state problem, massless states appear 
automatically because of the decoupling of the U(l) and SU(iV) degrees of freedom that 
constitute U(iV). More explicitly, we may introduce the U(l) operators 

a(k+) = ^tr[a(A; + )] and (3{k + ) = ^-tr[b(k + )}, (4.8) 
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which allow us to decompose any U(iV) operator into a sum of U(l) and SU(iV) operators: 

a(k + ) = a(k + ) ■ l NxN + a{k + ) and b{k + ) = (3(k + ) ■ l NxN + b{k + ), (4.9) 

where a(k + ) and b(k + ) are traceless N x N matrices. If we now substitute the operators 
above into the expression for the supercharge ( |4.7| ), we find that all terms involving the 
U(l) factors a(k + ), )3(k + ) vanish - only the SU(iV) operators a(k + ), b(k + ) survive, i.e. 
starting with the definition of the U(iV) supercharge, we end up with the definition of the 
SU(iV) supercharge. In addition, the (anti) commutation relations [aij(ki), a'fa)] = 
and {bij{ki) , ft {k 2 )} = imply that this supercharge acts only on the SU(iV) creation 
operators of a fock state - the U(l) creation operators only introduce degeneracies in the 
SU(iV) spectrum. Clearly, since Q~ has no U(l) contribution, any fock state made up 
of only U(l) creation operators must have zero mass. The non-trivial problem here is to 
determine whether there are massless states for the SU(iV) sector. We will address this 
topic next. 

4.2 The Proof for (1,1) Model. 

It was pointed out in the previous subsection that a zero mass eigenstate is annihilated 
by the light-cone supercharge (|4.7| ): 

Q-\^} = (4.10) 

We wish to show that if such an SU(iV) eigenstate is normalizable, then it must involve 
a superposition of an infinite number of Fock states. The basic strategy is quite simple: 
normalizability will impose certain conditions on the light-cone wave functions as one or 
several momentum variables vanish. Moreover, if we assume a given eigenstate |\l/) has 
at most n partons, then the terms in Q~\^>) consisting of n+ 1 partons must sum to zero, 
providing relations between the n parton wave functions only. We then show these wave 
functions must all vanish by studying various zero momentum limits of these relations. 
Interestingly, the utility of studying light-cone wave functions at small momenta also 
appears in the context of light-front QCD 3+1 M. 

In order to proceed with a systematic presentation of the proof, we start by considering 
the large N limit case. This simply means that we consider Fock states that are made 
from a single trace of a product of boson or fermion creation operators acting on the 
light-cone Fock vacuum |0). Multiple trace states correspond to 1/N corrections to the 
theory, and are therefore ignored. In this limit, a general state |\1/) is a superposition of 
Fock states of any length, and may be written in the form 

„=2r=0 P J ° Vqi-'-Qn 

f^ r \ qi ,..., qn )tr[c\ qi )...c\ qn )]\0), (4.11) 
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where c j[ (q + ) represents either a boson or fermion creation operator carrying light-cone 
momentum q + , and fp ,r denotes the wave function of an n parton Fock state containing 
r fermions in a particular arrangement P. It is implied that we sum over all such 
arrangements, which may not necessarily be distinct with respect to cyclic symmetry of 
the trace. 

At this point, we simply remark that normalizability of a general state |\1/) above 
implies 

P+ ^ • • • dqn S( qi + ■ • - + q n - P + )\f P n ' r \qi, . . . , g »)| a < oo (4.12) 

o gi • • • q n 

for any particular wave function fp' r . Therefore, any wave function vanishes if one or 
several of its momenta are made to vanish. 

We are now ready to carry out the details of the proof. But first a little notation. 
We will write \^( n ,m)) to denote a superposition of all Fock states - as in ( |4.11| ) - with 



precisely n partons, m of which are fermions. Such a Fock expansion involves only the 
wave functions fp , and the number of them is enumerated by the index P. For the 
special case \^( n ,o)) (i-e. no fermions), there is only one wave function, which we denote 
by /( n '°) for brevity: 

= f P+ d ^-dQn s{gi + +gn _ p+) /K0)(gi ... ftl ) te [ a t( fi )...ot(g B )]|0). (4.13) 
Jo ^Jq\---q n 

There is another special case we wish to consider; namely, the state |^/( n ,2)) consisting 
of n parton Fock states with precisely two fermions. If we place one of the fermions at 
the beginning of the trace, then there are n — 1 ways of positioning the second fermion, 
yielding n — 1 possible wave functions. We will enumerate such wave functions by the 
subscript index k, as in fjf 1 ' , where k — 2, 3, . . . , n. The subscript k denotes the location 
of the second fermion. Explicitly, we have 

l*(n, 2 )> = t l P+ %=^%! + " " " + .n ~ P + )X 
fc=2 ° y/Ql---Qn 

ft' 2 \q u ...,q k ,..., q n Mb\ qi )a\q 2 ) . . . b\q k ) . . . a\q n )}\0). 

(4.14) 

Of course, depending upon the symmetry, the n — 1 Fock states enumerated in this way 
need not be distinct with respect to the cyclic properties of the trace. This provides us 
with additional relations between wave functions - a fact we will make use of later on. 

Now let us assume that |\1/) is a normalizable SU(iV) zero mass eigenstate with at 
most n partons. Glancing at the form of (f4.7|), we see that the n + 1 parton Fock states 
containing a single fermion in each of the combinations Q~|\l/( ra) o)) and Q~\^( n ,2)) must 
cancel each other to guarantee a massless eigenstate. This immediately gives rise to the 
following wave function relation: 

gi+2g 2 ( 0) , x gi + 2g n+1 ( 0) 

; J [Qi + Q2,q3,- ■ -,qn+i) ; J [Qn+1 + 9l> 92, • • • , Qn) — 

9i + <?2 9i + q n +i 
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_oV (?1 V^ Qk+l-Qk An,2), . s. 

~ n h (qk + i + q k ) 3 / 2h (gi ' ■ ■ ■ ' qk ~ U ^ ft+1 ' gfc+2 ' " " gn+l) ' 

(4.15) 

In the limit g^ — > 0, for 3 < z < n, this last equation is reduced to 

-fi ' (9l> • • • > 9i-l> 9i+l; • • • ; 9n+l) 



1 ,(n,2), 



V9i+1 ' 
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/i_i (gi,...,ft_i,g»+i,...,g n+ i) = 0. (4.16) 
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An immediate consequence is that any wave function f\ n ' ' for i = 3,4, . . . ,n, may be 
expressed in terms of / 2 • Explicitly, we have 

fi U,2) (<ll,<l2,---,qn) = J — ft' 2 \<lli<l2, ■■■,9n), l = 3, 4, . . . , 71. (4.17) 

Moreover, the limit g 2 ~ * of equation Q4.15| ) yields the further relation after a suitable 
change of variables: 

/ (n ' 0) (9i,ft,g3,...,9n) = -Mf^fa, ft, 93, • • • , ft). (4.18) 

n V ^ 2 

Finally, because of the cyclic properties of the trace, there is an additional relation 
between wave functions: 

fi (Ql, 92, ■ ■ ■ , ft, • • • , Qn) = -fn-i+2(qi, ft+1, ■ ■ ■ , Qn, 9l, ft, • • • , ft-l)- (4-19) 

Setting i = 2 in the above equation, and z = n in equation (|4.17|) , we deduce 



/ 2 (n ' 2) (ft, ft, • • • , ft) = -J—f2 n ' 2) (q2, 93, • • • , ft, ft)- (4.20) 
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Combining this with equation ( |4.18| ), we conclude (— + : ^-)f^ n ' \qi, . . . , q n ) — 0, where 



we use the fact that the wave functions f( n >°) are cyclically symmetric. Thus f( n < ^ must 
vanish. It immediately follows that f^ vanish for all i as well. 

To summarize, we have shown that if |\1/) is a normalizable zero mass eigenstate, 
where each Fock state in its Fock state expansion has no more than n partons, the 
contributions |^( n ,o)) and \^( n ,2)) i n this Fock state expansion must vanish. Since we 
may assume |\1/) is bosonic, the only other contributions involve Fock states with an even 
number of fermions: |\I/(n,4)), |^(n,6)}, and so on. We claim that all such contributions 
vanish. To see this, first observe that the n + 1 parton Fock states with three fermions in 
the combinations Q~\^( n ,2)) and Q~\^( n ,4)) must cancel each other, in order to guarantee 
a zero eigenstate mass. But our previous analysis demonstrated that \^( n ,2)) = 0, and 
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thus the n + 1 parton Fock states with three fermions in Q~\^( n ,4)) alone must sum to 
zero. 

We are now ready to perform an induction procedure. Namely, we assume that for 
some positive integer k the state \^f ( n ,2[k-i])) vanishes. Then the n+1 parton Fock states 
in Q~\^f) which contain 2k — 1 fermions receive contributions only from Q~\^( n ,2k)) in 
which a fermion is replaced by two bosons. This has to sum to zero. We therefore obtain 
a relation among the wave functions fp by considering the action of the supercharge 
Q4.7D in which a fermion is replaced by two bosons. Keeping in mind that we are free to 
renormalize any wave function by a constant, we end up with the following relation: 

2 t fp ' ( s i) • • • ' s «-i' s i + s «+i' s «+2) • • • -i s n+i)- ( : rrj^ = 0. (4-21) 

p \ s i+l + s i) 

It is now an easy task to show that the wave functions fp appearing in equation ( |4.21|) 
must vanish; one simply considers various limits Sj — > as we did before. This completes 
our proof by induction. Namely, there can be no non-trivial normalizable massless state 
with an upper limit on the number of allowed partons. Of course, this proof is valid only 
in the large N limit. We now turn our attention to the finite N case. 

Let us define Q~[ ead to be that part of the supercharge Q~ that replaces a fermion 
with two bosons, or replaces a boson with a boson and fermion pair. As in the large N 
case we begin by assuming that we have a normalizable zero mass eigenstate |\1/) which 
is a sum of Fock states that have at most n partons. The proof for finite N consists of 
two parts. First, we consider bosonic states consisting of only n parton Fock states that 
have at most two fermions, and show the wave functions must vanish. We then invoke 
an induction argument to consider n parton wave functions involving an even number of 
fermions, and show they must vanish as well. 

The additional complication introduced by the assumption that N is finite is that a 
given Fock state may involve more than just a single trace. However, note that Q~[ ead 
cannot decrease the number of traces; it can either increase the number of traces by 
one, or leave the number unchanged. Thus we have a natural induction procedure in the 
number of traces as well. Since the terms in Q~[ ead have only one annihilation operator, 
it acts on a given product of traces according to the Leibniz rule: 

Qg ad (trL4]tr[£] . . .) |0> = (Qj^tr[A]) tr[B] . . . |0> + 

{-l) F{A hr[A]Q Tead {tr[B]...m. (4.22) 

Schematically, the general structure of an arbitrary Fock state with k traces has the form 
yKMiA,..,i*) tr [( 6 t)ii t . . . t] . . . tr [(fct) <* a t . . . a t]| ), (4.23) 

where n denotes the total number of partons in each Fock state, and the integers ii,i 2 , . . . 
denote the number of fermions in the first trace, second trace, and so on. We will always 
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order the traces so that the number of fermions in each trace decreases to the right. The 
index P labels a particular arrangement of fermions. 

We now consider the n + 1 parton Fock states of Qi~ ea d\^) that have precisely one 
fermion. The only possible contributions involve three types of wave functions; f^ n,Q \ 
fp and /( n ' 1 ' 1 ) (we only include the permutation index P if there is more than one 
distinct arrangement). If these three wave functions contribute to the same one fermion 
Fock state, then the distribution of bosons in the Fock state corresponding to fp 
determines the distribution of bosons for /( n >°) and /w 1 ' 1 ). We allow Q~[ ead to act only 
on the first trace in both /( n >°) and fp , and only on the second one in f\ n > 1 ' 1 ). If there 
are more than two traces in these states they must be identical in all the components, 
and so don't play a role in the calculation. Thus, it is sufficient to consider states with 
two traces only. Such a state has the form 
f d m+n q 



£(n+m,0) 

xtr 



(9ii 

a\q x ) 



y/Ql ■ ■ ■ Qn+m 
■ ■ ■ i Qm\Qm-\-li 



%i H hq. 



n+m 



-I* 



j Qm+nj 



a\q n 



tr 



a\q m+ i)...a ! (q m+n ) |0) 



rP + d m+n ~ 2 qd Pl dp 2 .. 
-I <?(gi + 



10 y/qi-.. q U A 

f(n + m,l,l) {puqu _ 



m-2PlP2 

■ , Qm\P2, Qm+3, ■ ■ ■ , Qm+n) x 



+ q n+m - 2 +P1+P2- P" 



xtr 



b\pi)a ] (qi) . . . a\q m )]ti{b ] (p 2 )a) {q m+ z) . . . a\q.. 



m+n, 



+ 



+ J2jP m ' (PU P[qi ■ ■ ■ qm-2] P2]\q m +1 ■ --qm+n) X 
P 

xtr (b\pi)P[a ] (qi) . . . a f (g m _ 2 ); 6 f (p2)]) tr [a f (g m+ i) . . . a ] (q m+n )\ \ |0) 



(4.24) 



where we have summed over the index P representing all possible permutation arrange- 
ments of bosons and fermions that contribute. We then find: 



F(p, q 1 






qm+2 qm+1 x(n+m,l,l) 



(q m +2 + qm+i) 3/2 ' 



qm+n ) 1 
[P, q\ . . . q m \qm+l + qm+2 j qm+3 



(4.25) 



qm+n) 



where F is the contribution from f( n+m '°) an d fp m ' . Now we see that the limit 



q-m+i —> gives: f( n + m ' 1 ' 1 ) = 0. Thus if Q4.24J) represents a contribution to the massless 
eigenstate state |^/), then |$) takes the form 

rP + d m+n-2 qdK + 



m 



y/qi ■ ■ ■ q n + m -2 



S(qi H \-q 



n+m— 2 



(P + ~ K+)) 



36 



f p+ dq m _ 1 dq m . 

•'O ^qm-lQm 

f {n+m,0) (qi, ■■■, q m \q m +i, • • • , 9m+n)tr(a f (gi) . . . a f (g m )) 
[ p+ dpidpi + 

"'o VP1P2 

2 /p n+m ' 2) (^l> P [?l> • • • » Qm-2] Pi] km+1, • • • , gm+n) (4.26) 

P 

tr(6t(p 1 )P[o t (gi) • • • a^-a); 6 t (p 2 )])tr(at(g m+1 ) . . . a\q m+n ))] |0> 

and Q z ~ arf acts only on the terms in the square brackets. All these terms have only one 
trace, which is a scenario we already encountered in the large N limit case. Using the 
results of that discussion, we find that the only massless solution of the form Q4.26 ) is 



the trivial one. This is the starting point of the induction procedure for finite N. 

As explained earlier, we look for n parton Fock states in the expansion for |\P) that 
have 2k fermions (k > 1), To finish the proof we need to show that for any k the only 
allowed wave function is the trivial one. ^From the large N result we know there are no 
such one trace states. We now consider the state with an arbitrary number of traces, 

r P+ A a Ac 

l*W*)> = £/ ■ dSn 6(s 1 + --- + s n -P+) (4.27) 

p JO y/Si ■ ■ ■ S n 

fp 2k \ Sl ...s ix \...\... Sn )tr ( c t( Sl ) . . . c t( flil )) tr (. . .) tr (. . . c t(a B )) |0), 
then the analog of ( f4.21|) for such states reads: 



22 fPi' ( S l- ■ ■ \ S ja ■ ■ ■ S i-1> S i + S i+l5 S i+2 ■ ■ ■ s ja+k a \ ■ ■ ■ S n+l)- ( j r^ — 0. (4.28) 

i \ s i+l "T Si) I 

Here, J2'i rneans that for each trace we should include one additional term with "z" = 
3a + k a , "i + 1" = j a if c corresponding to both j a + k a and j a is a. If the number of 
traces is a, we introduce 

a-l 

Ja = J2 k b- 
6=1 

If any of the blocks tr(. . .) in the state for which (|4.28|) is written contains two or 
more fermions, then, as in the large N case, all the corresponding wave functions fp 
vanish. So we only need to consider the states of the form: 

|*(n,fci+i,...)) =J2 / d P d qfp M+1, '" ) (P^q^---^k 1 \P2,qk 1 +i,---,qk 1 +k 2 \---) x 
p J 

tr {b ] (p x )a ] (q x ) . . . a f (g fcl )) tr {b ] (p 2 )a ] (q kl+1 ) . . . a ] (q kl+k2 )^j . . . |0). 

(4.29) 
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Let Q denote that part of the supercharge Q~ which replaces a fermion with two bosons. 
Let us consider the result of such a change in the first trace. Suppose there are a traces 
having the same form as the first trace. Then without loss of generality, we may assume 
they are the first a traces. Then using the symmetries of the wave functions we find: 

l ^ r p+ 



| r-f ' 

0|*(n,fci+i,...)) = -^7k=Y1 / dkdpdq 
ly In P Jo 



p 

<.(n,fci+l,...) / i i \ v^ Pi ~ * J- / i\6+l 

fp {pi,Qi,---,Qk 1 \P2,qk 1 +i,---,g2k 1 \---)2^ / (~ X ) x 

6=i Pb ^k(p b - k) 



tr [b^ (jh)^ (qi) ■ ■ ■ a ] (q kl )j . . . tr (a t (k)a f (p b - k)a ] (q {b _ 1)kl+1 ) . . . a\q bkl )j . . . |0) 

1 f p+ p-i — 2k 1 / \ 

Y^ / dkdpdq =tr (aUk)aUpi — k)aUqi) . . . aHgt,)) x 

d JO V\ Jhfr,-, - h\ V ' 



2^/2it^Jo ' ' pi ^Jk(p 1 - k) 

tr (b\p 2 )a\q kl+l ) . . . a\q kl+k2 )) . . . |0> E(-l) 6+1 (-l) 6+1 x 



r(n,ki+l,...) , I I \ 

Jp [Pi,Qi, ■ ■ ■ , Qk!\P2, qk!+i, ■ ■ ■ , Qkx+kil ■ ■ ■)■ 

If the above expression vanishes then the only solution is the trivial one in which all wave 
functions vanish. This finishes the proof of the induction procedure for the finite N case. 
The extension of the proof to massive bound states is straightforward. Firstly, assume 
|\&) is a normalizable eigenstate of 2P + P~ with mass squared M 2 ^ 0. Then, since 
p- = -j^(Q-) 2 , the state 

|*) = |tt) + a Q-|tt) (4.30) 

for a 2 = \/2P + jM 2 is a normalizable eigenstate of the supercharge Q~ , with eigenvalue 
\Jol. We therefore study the eigen-problem Q~\ty) = -\^f). The resulting constraints on 
the wave functions may be obtained by modifying our original expressions by including a 
wave function multiplied by a finite constant. However, in our analysis, we always need 
to let some of the momenta vanish, and therefore this additional contribution vanishes. 
The analysis (and therefore the conclusions) remains unchanged. 

We therefore conclude that any normalizable SU(iV) bound state (massless or massive) 
that exists in the model must be a superposition of an infinite number of Fock states. 

4.3 Higher Dimensional Theories. 

In this subsection we extend our theorem to the two-dimensional supersymmetric theories 
obtained as the result of dimensional reduction from D > 3 dimensions. The most 
important cases are D — 4, 6 and 10 which have 2, 4 and 8 supersymmetries in two 
dimensions. Below we consider only large N case, the generalization to arbitrary SU(N) 
group is trivial repetition of the arguments given in previous subsection. 
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Again our starting point is the fact that if there is normalizable eigenstate of Hamil- 
tonian having finite length than its main symbol satisfies the condition: 

QTeJy >= 0, (4.31) 

where Qi~ ea d * s the part of supercharge increasing the number of partons. In three dimen- 
sional case we had only one supercharge Q~ , for general D dimensional SYM reduced in 
1 + 1 there are D — 2 supercharges, each of them squared gives P~ and ( |4.31| ) should 



be true for all of them. In general different supercharges are not anticommute with each 
other, but since we consider quantization near trivial classical configuration (with no 
monopoles and no external charges) then they do. It is easy to derive the general form 
of supercharge: 

QZ = I ™(b$(k)J ij (-k)-(J ij (-k))%(k))+ (4.32) 



+ ^= IJkM^A^AjUk^i-k), 
J ij(~ k ) = 7T1k= d P i , \ ( a ki(pWkj( k + P)~ ajl(p)alk( k + P)) + ( 4 - 33 ) 



\ fk b — Or, 1 fk 

+ — = dp n ==J=al k {p)a{ j (k - p) + ^= dpbUp)bU k - P) + 
2^2-nJo yfp(k-p) V2tt Jo 

+ 4= f dp (b a J(p)b a kj (k + p)- bf k {p)bUk + P 



. 2tt 

In the above expression we introduced d = D — 2 kinds of bosons (I = 1, ...,d) and d 
kinds of fermions (a = l,...,d) which we get as the result of compactification. The \x 
is nonzero constant depending on D and M are combinations of d dimensional Dirac 
matrices: 

Mfj = (7/tJ - Ul T i U- (4-34) 

As before our proof is based on the induction on the number of fermionic operators 
in the state. First we consider main symbol being superposition of purely bosonic states 
and ones containing two fermionic operators. Now we have d types of bosons and d types 
of fermions so some additional indices should be included in the wavef unctions. Defining 
bosonic indexes to be capital letters A, B... and fermionic ones to be Greek letters we 
write: 

n P+ %%to + .- +dn- P + ) E /[1.^1-Sn) X 

o y/qi...q n A 

x tr[o^(g 1 )...oi n (g„)]|0>, (4.35) 

k=l JU yHl---Hn A }0i 

x ^[4 1 (9i)...4 fe „ 1 fe-i)&i 1 (%)al 4fe fe+i)--- a A n _ 2 (9n-i)fei 2 (g„)]|0> . (4.36) 
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It is now easy to find the one fermionic part of the result of action by ( |4.32| ) on the 



main symbol of the state. The vanishing of this contribution leads to the generalization 
of the equation ( (4.15 ) : 



n ( 2q n +p (0) 2qi±PAo) (n , nn „ Ol 

Oa/3- : /Ai...A„Wl--9n-l,9n + P . Ji 1 ...A n 9l + P—Qn-1, Qn) ~ 

P \ Qn+P qi+p J 

2 ^ qk+i — qk j- f (2)k , . \ . 

■J= 2^ ( + qk y/2 d A k A k+1 J [Al ...A k _ 1 aA k+2 ...A n p]{ ( ll-- ( lk-l,(lk + 9fc+l, 9fc+2-9n,Pj + 

\ qn ~r p yi "t f / 



2 /^ y^ M T k+1 A k ,(2)fe 



E Afc+ ^ fc / [ A ) 1 fc .A fc _ l7 A fc+2 ...A n /3](gi-gfc-i^ gfc + q k +i,qk+2-qn,p) = (4.37) 



This equation should be true for any possible Ai...A N , a and /3. We will show that the 
only solution of such system of equations is trivial one so all the /> j and /r j vanish. 
This will be proven by induction. First we note that if Ax — ... — A n and a = /3 then 
equation Q4.37Q is reduced to ( |4.15| ) written for u A a , and ft A AaA A , and as we saw 



this leads to 

f[A...A] = °> f[A...AaA...Aa] = ° ( 4 - 38 ) 

for arbitrary A and a. The next case to consider is Ax = ... = A n , a ^ (3. Using relation 
just found the ( |4.37| ) for this case again gives us ( [4.15| ), but this time correspondence 



reads: 



/^-J^/gL 



.A] 

'[A. ..AaA.. .A/3] 

The proven property of fl4.15|) together with trivial identity 



f (2)k - fll.A aA ...A8V (4-39) 



£ M p CA Mf B = 4dS BC (l - 5ac) (4.40) 

a/3 

leads to 

f[BA...A] = 0' f[A...AaA...A/3] = ® (4-41) 

for any A, B, a, (3 (A could be equal to B and a to 0). We use this equation as starting 
point of the induction procedure. 

Let us introduce one useful function. For each set {Ai...A k } we define n({Ax...A k }) 
to be the maximal number of identical A in the set: 

n{{Ax...A k }) = max fc 6(A t - I) j . (4.42) 
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/£Lw = ° «/ n({A-..A n })>n-l 



In terms of this new function our result ( 4.41 ) can be rewritten as 

,(0 
'[A 

f{Z k ..A h _ ia A k ...A n . 2 p] = if n{{A l ...A n }) = n-2. (4.43) 

This condition will be used as starting point of induction then the assumption of induction 
procedure is: 

/[2L4n]=° «/ n({A l ...A n })>m 

/SU-M^,fl = (1 </ n({A...AJ)>m-l (4.44) 

and we checked it for m — n — 1. In our induction procedure we will decrease parameter 
m instead of increasing it. To perform the proof, we start with writing Q4.37 ) for the set 
{y4i...y4fc} with n = m: 

2 ^-v qk+i — Ik x ,(2)fc / . \ _ 

nil \ ^f-ftl An - xB {<li--<ln-iAn +p) ~ ^^fSL.AAQi+P-Qn-uqn)) ■ (4.45) 
\q n + P Qi+P I 

If a — 3 then the right hand side is zero and we have a recurrent relations for f,J A nA . . 
with different k. Due to the presence of 5 symbol these relations would connect only f^ k 
inside some clusters (for m < n) and the boundary elements of such clusters should be 
zero. Thus we deduce that if n({Ai...A n ^ 2 }) = rn — 2 then 

flZ.A k -, a A k ...A^ a] = 0- (4-46) 

For a t^ /3 we consider different limits qi — > in ( |4 . -1 -5| ) : 

~ V -SA 1 A 2 fiaA 3 ... An p ] (q2~-qn,p) = -n[iMf lB ff M „. An (p...q n -u q n ) (4.47) 



1 ? (2)i 



jqi 

for % — 1 and 



^A i A i+ i/[ J 4 1 ... J 4 i _ 1 ai4i +3 ...>l n ^]Wl—9*-l'5i+l— 9n,P) 



I - e {2)i 



/ ^S Al A^ 1 f[A 1 l ...A l . 2a A 1+1 ...A n l3]( ( ll- ( li-U^+l- ( ln,p) (4.48) 

y/qi-l 

for 1 < i < n. Since n < n — 1 then there exist z < n: Saa- 1 = then from above 
equations we deduce for a ^ (3: 

f[Z..A n] = 0, ^({A-A}) = m - 1 
f[Z..A k ^A k ...A n . 2P] = 0, ^({A-AJ) = m - 2. (4.49) 
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This finishes the proof by induction. Thus we have proven that equation ( |4.37| ) doesn't 
have any normalizable solutions. 

To show that there are no finite length bound states we now turn to the analog of 



equation (|4.21| ). This analog reads: 



7 , Aifp. (Si...Sj_i, Si + Si+i, Si+2-"Sn+l, 



x 



x ( s Xk +1 ( Sl+ \ w 2 + » M %i w r-^- ) = o, (4-50) 

\ I+1 (S i+1 + Si) 6 / 2 l+1 y/S i+1 + Si J 

where P, describes different permutations of A and a. This equation gives linear relations 
between wavefunctions inside the block of a in 

tr(...6 t a t ...a t 6 t ...)|0 >, (4.51) 

"boundary" elements (when index i corresponds to fermions) vanish, so as in the case 
3 — > 2 all the /( 2fc ) are zero. This completes the proof for general compactification. 

4.4 Bound States in DLCQ. 

In the previous subsection we proved that the continuum formulation of the theory does 
not have any normalizable bound states with a finite number of partons. Our proof 
used the behavior of wave functions at small momenta arising from the normalizability 
assumption. Neither of these properties can be used in DLCQ, however. Here we consider 
some simple examples of massless DLCQ solutions with n bosons to help shed some light 
on the relation between DLCQ solutions and the solutions of the continuum theory. For 
simplicity, we work in the large N limit case. 

We write the momentum of a state in DLCQ in terms of the momentum fraction qi 
where qi = —P + , and the r^ are positive integers. The wave function of such a state 
is /^ n '°- ) (r 1 , . . . ,r n ). There are two conditions that must be satisfied to show that it is 
massless. One is the that the coefficient of the term with one additional fermion that is 
produced by the action of Q~ is zero. This condition gives the relation, 

*±±lfM( ri , . . .,r n . 1 ,r n + t)- *"=!±?/M) (ri , . . . ^ x + t,r n ) = 0. (4.52) 
r n + t r n _i + 1 

where t correspond to the momentum fraction of the one fermion. The second is that the 
coefficient of the state with two fewer bosons and one additional fermion which is also 
produced by the action of Q~ is zero. This condition gives the relation, 

t — 2k 

£ uT^A f{nfl) fa' • • • ' r «" 2 ' fc ' * - Wm-i-nw) = 0. (4-53) 

For the case where all r, = 1, and the total harmonic resolution is n, it is trivial 
that eqn( [4.52] ) is satisfied since there is not enough resolution to increase the number of 
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particles in the state. It is also easy to see from eqn ( |4.53| ) since the coefficient of the 
one term in the sum is zero. Thus the wave function /( n ' )(l, 1, ....1) is a massless state 
for every resolution 

To discuss additional solutions it is useful to start by considering eqn( [4.52|) . The case 
t = 1, gives the equation 

/^°)(r 1 ,...,r n _ 2 ,r n _ 1 ,r w + l) = 2 ;- 1 + 1 Tn ± ^ - /<"»°>(n, ■ ■ ■ , r n ^r n _, + 1, r n ). 

LT n + 1 T"n—1 * t 

(4.54) 
This equation is trivial to satisfy if r $ = 1 for all i. The contributions in eqn( [4.53| ) come 
from the two terms in the sum, k = 1, t = 3 and k = 2, t = 3. Each term has the 
same coefficient but of opposite sign and cancel. Therefore the state /( n, °)(l, ...1,2) is a 
massless state for all resolutions 



The next case t = 2 in eqn(|4.52|) gives 



/^°)(r 1 ,...,r w _ 2 ,r„_ 1 ,r w + 2) = 2 ^- 1 + 2 r " + 2 - /<"»°>(n, ■ ■ ■ , r n . 2 ,r n ^ -f 2. r 



i ' nj 



2r n + 2 r n _ x + 2* 

(4.55) 
Using ( 4.54 ) twice we find: 



/^)(r 1 ...r w _ 2 ,r w _ 1 ,r„. + 2) = 2rn ~ 1 + 1 r J^ f(^)(r 1 ...r^,r n . 1 + l,r n + l) 

Lr n + o v n _i + l 

2r n _! + 1 r n + 2 2r n _i + 3 r n + 1 c{nfi) 



2r n + 3 r n _! + l 2r n + 1 r n _! + 2 
Comparing with ( [4.55 ) we have: 



/^(r 1 ...r n _ 2 ,r n _ 1 + 2,r n ). (4.56) 



/ Cn.-.r^r^ + ^r^l + + I U. (4.57J 



Using relation (|4.52|) several times we can always express an arbitrary wave function in 



the following form: 

/^'°Vi...r n ) = C( ri ...r n )f^°\l...l, L + 1,1) (4.58) 

where L = r\ + ... + r n — n and C(ri...r n ) is some nonzero coefficient. The two massless 
states we found above correspond to L = and L = 1. Choosing r\ = ... = r n _ 2 = r n = 1 
in ( |4.57| ) we find, 

/ (n ' 0) (l...l,(L-l) + 2,l) =0 for L>2 (4.59) 

due to monotonic behavior of the function in the parenthesis. Then using (|4.58|) we 
conclude that all the wave functions with L > 2 vanish. So the only case we need 
consider is L = 2. In this case ( |4.52| ) has only two nontrivial cases: t — 1 and t = 2 
which are given by (|4.54j) and Q4.55|) . In the second of these equations we can only have 
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r\ — . . . — r n — 1 so it is trivially satisfied. Equation ( |4.54| ) however gives a nontrivial 
relation for the wave function: 

/M)(l,...,l,2,2)=/^°)(1,...,2,1,2) = ...= 

/^°)(2,...,1,1,2) = ^°)(1,...,1,3). (4.60) 

finally we must show that eqn( |4.53| ) is satisfied which is straight forward. 

These are only a few examples of massless states, and there are in fact many more 
in DLCQ ||. But the results of our numerical analysis show that the states we just 
described are closely connected with the massless states in the continuum. Let us for- 
mulate this relation for N = oo. In this case only single trace states should be kept in 
the spectrum and DLCQ massless states have the following structure. The state first 
appear as tr(a^(l) . . . a^(l)) at resolution P, then one can trace it to resolutions P + 1 
and P + 2 as states with wavefunctions ( |4.59| ) and ( f4.60| ). As we just proved at higher 



resolutions there are no massless states containing exactly K partons, however at any 
resolution K > P there is exactly one massless state whose wavefunction is localized pre- 
dominantly in the sector with P partons (V . So it is natural to collect all such states in 
the single sequence and to call the limits of this sequence " continuum massless state with 
P bosons" , although as we saw the wavefunction of continuum state has contributions 
from sectors with different number of partons. The interesting feature of this theory is 
that such "continuum massless states with P bosons" are the only bosonic massless states 
seen by DLCQ (in principle the theory in the continuum might have massless state whose 
wavefunction is localized in sector with infinite number of partons, but we will ignore 
this possibility). Thus one can easily count bosonic massless states at any resolution P: 
they are just images of states with P bosons for all P < K, thus there are K — 1 such 
states. Acting on any of such states by Q + we can get the fermionic massless state (then 
there are also K — 1 of them), while acting by Q~ doesn't give any new state (the result 
is zero). We will do the counting for finite A" case in the next subsection. 

In the continuum limit we have proven that there are no massless normalizable states 
with a finite number of particles. However, at each finite value of the harmonic resolution, 
one obtains an exactly massless bound state, but as the harmonic resolution is sent to 
infinity, the number of Fock states required to keep the bound state massless must also 
be infinite. 

4.5 Counting of Massless States in DLCQ. 

Finally in this section we will count the massless states in DLCQ as function of resolution, 
keeping the number of colors A" finite. However we will assume that A" is not too small, 
so that the relation between A" and resolution K: K < N 2 — 1 is satisfied. We will need 
this condition in order to insure all states of the form tr(c^ . . . c 1 ) . . . tr(c' . . . c^)|0) are 
linearly independent unless they are related by either cyclic permutations in one of the 
traces or permutations of traces themselves. The simplest example of violation of this 
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condition is the state 

tr(a t (l)a t (l)a t (l))|0)=0 

for SU(2) (here 3 = 4—1). Although for K > N 2 — 1 some conclusions can be made, 
the different N and K requires special consideration and we are not going to proceed in 
this direction. From the numerical perspective we should mention that in our calculation 
K < 11, so the only excluded values of N are 2 and 3. 

As soon as the condition K < N 2 — 1 is satisfied the DLCQ Fock spaces for SU(N) 
and SU(oo) are the same if all multitrace states are taken into account. Moreover our 
numerical analysis strongly suggests that the number of massless states is the same for all 
iV > \/K + 1, while wavefunctions depend on N. However talking about SU(oo) Fock 
space one usually considers only single trace states as fundamental ones, while multitrace 
states are though of as the system of free bound states. Let us explain the reason for 
this. The light-cone Hamiltonian can be written in the following schematic form: 

1 11 1 

p" = — P~ = acjjCij + —pcljiccc)^ + jjPidc^ijicclij + — ^(cW^c)^. (4.61) 

The jj is introduced in order to make the eigenvalues of p~ finite as N — > oo. Let us 
consider two eigenstates of p~ , which are chosen to be combination of single traces in the 
large N limit: 

p~A\0) =m A A\0), (4.62) 

p-B\0) =m A B\0). (4.63) 

This is equivalent to the following commutation relations: 



[p ,A}= m A A + — Y^ /Mntr(c f . . . c j c) + — ^ v A ti(c ] . . . c ] cc) + O ( — 
p~,B = m B B + — J2 Atijntr(c t . . . c ] c) + — ^ z/ B ntr(c f . . . c ] cc) + O ( -^ 



(4.64) 



We introduced a convenient notation for the normalized trace here: 

ntr( c t ...c t c...c ) = ^^(c f . . . c%(c.. . c) y . (4.65) 

n 

This way the state ntr(c^ . . . c')|0) has a finite norm in the large A^ limit. ^From the 
equations (f4.64|) one can easily see that 



p~AB\0) = {m A + m B )AB\0) + O (J^ , (4.66) 

i.e. we indeed have a combination of two free states in the large iV limit. In particular 
this fact may be applied toward the classification of DLCQ massless states at iV = oo: 
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the multitrace state is massless if and only if all of the traces involved correspond to 
massless states. We also mention the trivial fact that if state A has resolution Ka and 
B has Kb then AB is the state at resolution Ka + Kb- 

Let us summarize what we have learned so far. The number of massless states in 
SU(N) theory at resolution K < N 2 — 1 is the same as one for SU(oo) theory if the 
multitrace states are included in the latter. On the other hand due to the fact that 
multitrace massless states in SU(oc) have special structure (namely any single trace in 
them is massless state itself), their number can be calculated from the known number of 
massless states written as linear combination of single traces. The remaining part of this 
subsection is devoted to such calculation. 

As we found in the end of the last subsection there are 2(K — 1) single trace massless 
states at resolution K. We will show how this information can be used in order to count 
the total number of massless states. Let us introduce some notation first. The value we 
want to calculate is iV* — the number of massless states at resolution k. We also define 
Nl m as number of such massless states at resolution k that the resolution of any single 
traces in them is greater or equal to m. Then for example N^ = Nf, and N^ = 2{k — 1). 
Finally we define / n (m) to be the number of different massless states containing n traces, 
each of which corresponds has resolution equal to m. Then one can derive the recurrent 
relation: 

[k/m] 
N ( m) = ^+1) + ^ Um) ^+1) + ^J (4 6?) 

n=l 

The starting point of the recurrent procedure are the relations N k = 2{k — 1) and 
N^ = for m > k. In order to apply ( |4.67|) we only have to evaluate f n (m). This will 
be our next task. 

To calculate f n (m) we first assume that we have only bosonic traces at our disposal. 
Let us count the states which contain p such traces. After combining identical traces 
together one can reduce the problem further by considering only states with special trace 
structure. Namely we will concentrate our attention on the massless states having rti 
traces of type i, for different values of i = 1 . . . r. Without the loss of generality we can 
require n% > - . . > n r > 1, one can also see that the relation J2 n i — P holds. If all rii 
are different then the number of massless states with fixed structure is given by simple 
formula: 

g(m) (g(m) - 1) . . . (g(m) - r + 1) , 

where g(m) = m — 1 is the number of bosonic single trace massless states. In general one 
gets additional combinatoric coefficient C{n\, . . . ,n r ) in the last expression, it is defined 
by the following rules: 

C(m,...,ni,n i+ i,...,n r ) = C(ni,...,ni)C(n i+1 ,...,n r ), if ^ > n i+1 , 

C(n,...,n) = 1 (4.68) 

" v ' a! 



4(3 



Now let us include fermionic traces in the picture. The only difference between them and 
bosons is the Pauli principle, so considering the product of q fermionic traces one has to 
choose all of them to be different. Thus the coefficient C for this case is 



C F (q) = C(l, 



1 



(4.69) 



Collecting all the information together we finally get: 

n i 

fn(m) = Y. —9(m) {g{m) - 1) . . . (g(m) - q + 1) x 
<j=o ?■ 

n—q 

x ^2 F{n — q, r)g(m) (g(m) — 1) . . . (g(m) — r + 1) 



(4.70) 



r=0 



F(p, r) 



J2 C(ni,...,n r ), 

{ni,...,n r } 



m + . . . + 7l r = p 
n\ > . . . > n r > 1. 



(4.71) 



Now we can use the relations ( 4.70Q , (|4.71 ) and (|4.68|) to determine all the coefficients 
f n (m) and then substituting them to the recurrent relation (|4.67|) one can find all the 



r (m) 



r(2) 



Nfc"'' . This in turn leads to the results for Nk — N% . Although we were not able to find 
an analytic expression for Nk as function of k, the number of states can be evaluated 
numerically for arbitrary k using the procedure we just described. We performed such 
calculations using Mathematica and the results for lowest resolutions are summarized in 
the Table(I]. For instance one can see that up to resolution 5 Nk = 2 k ~ 1 , but at higher 
resolutions this relations holds only approximately. 



k 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


N k 


2 


4 


8 


16 


32 


60 


114 


212 


384 


692 


1232 



Table 1: Number of multitrace massless states as function of resolution. 



5 Correlation Functions in SYM and DLCQ. 

The bound state problem we have studied so far is the traditional one for DLCQ. How- 
ever this is not the only calculation that can be done using this method. The problem 
of computing of correlation functions, more traditional for conventional quantum field 
theory, can also be addressed in the light cone quantization. Unlike the usual methods of 
QFT the results of DLCQ calculations are valid beyond the perturbation theory and thus 
they can be used for the testing the duality between the gauge theory and supergravity. 
There has been a great deal of excitement during this past year following the re- 
alization that certain field theories admit concrete realizations as a string theory on a 
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particular background |nj]. By now many examples of this type of correspondence for 
field theories in various dimensions with various field contents have been reported in the 
literature (for a comprehensive review and list of references, see 0). However, attempts 
to apply these correspondences to study the details of these theories have only met with 
limited success so far. The problem stems from the fact that our understanding of both 
sides of the correspondence is limited. On the field theory side, most of what we know 
comes from perturbation theory where we assume that the coupling is weak. On the string 
theory side, most of what we know comes from the supergravity approximation where 
the curvature is small. There are no known situations where both approximations are 
simultaneously valid. At the present time, comparisons between the dual gauge/string 
theories have been restricted to either qualitative issues or quantities constrained by 
symmetry. Any improvement in our understanding of field theories beyond perturbation 
theory or string theories beyond the supergravity approximation is therefore a welcome 
development. 

We will study the field theory/string theory correspondence motivated by considering 
the near-horizon decoupling limit of a Dl-brane in type IIB string theory |)5U|| . The gauge 



theory corresponding to this theory is the Yang-Mills theory in two dimensions with 16 
supercharges. Its SDLCQ formulation was recently reported in ||. This is probably the 
simplest known example of a field theory/string theory correspondence involving a field 
theory in two dimensions with a concrete Lagrangian formulation. 

A convenient quantity that can be computed on both sides of the correspondence 



is the correlation function of gauge invariant operators [|T| pq] . We will focus on two 
point functions of the stress-energy tensor. This turns out to be a very convenient 
quantity to compute for many reasons that we will explain along the way. Some aspects 
of this as it pertains to a consideration of black hole entropy was recently discussed in 



i^| . There are other physical quantities often reported in the literature. In the DLCQ 
literature, the spectrum of hadrons is often reported. This would be fine for theories in a 
confining phase. However, we expect the SYM in two dimension to flow to a non-trivial 
conformal fixed point in the infra-red |5^, |31]. The spectrum of states will therefore 
form a continuum and will be cumbersome to handle. On the string theory side, entropy 



density [p|]] and the quark anti-quark potential |24], [73], ^2| are frequently reported. The 



definition of entropy density requires that we place the field theory in a space-like box 
which seems incommensurate with the discretized light cone. Similarly, a static quark 
anti-quark configuration does not fit very well inside a discretized light-cone geometry. 
The correlation function of point-like operators do not suffer from these problems. We 
should mention that there exists interesting work on computing the QCD string tension 
T5| , |16j directly in the field theory. These authors find that the QCD string tension 



vanishes in the supersymmetric theories which is consistent with the power law quark 
anti-quark potential found on the supergravity side. This section is based on the results 
of paper M . 
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5.1 Correlation Functions from Supergravity 

Let us begin by reviewing the computation of the correlation function of stress energy 
tensors on the string theory side using the supergravity approximation. The computation 



is essentially a generalization of [11, 33 . The main conclusion on the supergravity side 



was reported recently in [42| but we will elaborate further on the details. The near 



horizon geometry of a Dl-brane in string frame takes the form 



,,/ I , U * dx\ + V 647r 9 ™ N dU 2 + Ju^g 2 YM NUdVt 2 



^9 2 ym r n \i MN Y ■ (5.1) 



'64tt 3 ' 

~[/ 6 / 



In order to compute the two point function, we need to know the action for the diagonal 
fluctuations around this background to the quadratic order. What we need is an analogue 
of [Q for this background which unfortunately is not currently available in the literature. 
Fortunately, some diagonal fluctuating degrees of freedom can be identified by following 



the early work on black hole absorption cross-subsections |56| [40fl . In particular, we can 
show that the fluctuations parameterized according to 

ds 2 = (l + f(x°, U) + g(x°, U)) g 00 (dx ) 2 + (l + 5f(x°, U) + g(x°, UJ) gil {dx 1 ) 2 

+ (l + f(x°, U) + g(x°, UJ) 9uu dU 2 + (l + f(x°, U) - ^g(x°, C/)) g w dtt 2 

e* = (l + 3/(x°,f/)-^(x°,f/))e 00 (5.2) 

will satisfy the equations of motion 

f"(U) + 7 -f { u)- Mw39 ™ Nk2 W) = o 

9"(U) + y(V) - ^g(U) - ^ g ™ Nk \ (U) = (5.3) 

by direct substitution into the equations of motion in 10 dimensions. We have assumed 
without loss of generality that these fluctuation vary only along the x° direction of the 
world volume coordinates like a plane wave e tkx . The fields f(U) and g(U) are scalars 
when the Dl-brane is viewed as a black hole in 9 dimensions; in fact there are the min- 
imal and the fixed scalars in this black hole geometry. In 10 dimensions, however, we 
see that they are really part of the gravitational fluctuation. We expect therefore that 
they are associated with the stress-energy tensor in the operator field correspondence of 
PH , PB| . In the case of the correspondence between J\f = 4 SYM and AdS 5 x S 5 , super- 



conformal invariance allowed the identification of operators and fields in short multiplets 
|JU| . For the Dl-brane, we do not have superconformal invariance and this technique is 
not applicable. In fact, we expect all fields of the theory consistent with the symmetry 
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of a given operator to mix. The large distance behavior should then be dominated by 
the contribution with the longest range. The field f(k°,U) appears to be the one with 
the longest range since it is the lightest field. 

The equation ( |5.3|) for f(U) can be solved explicitly in terms of the Bessel's function 



f(U) = U- S K 



3/2 1 



l^ 3 g 2 YM NU 



-2 



k). 



(5.4) 



By thinking of f(U) in direct analogy with the minimally coupled scalar as was done in 
EH |83|, we can compute the flux factor 



T 



lim 

U 



1 



oo 2«f 



y/99 



UU e -2( 



^dulog(f(U)) 



U=U 



NU 2 k 2 
^9ym 



N 3 / 2 k 3 
^9ym 



+ 



(5.5) 



up to a numerical coefficient of order one which we have suppressed. We see that the 
leading non- analytic (in k 2 ) contribution is due to the k 3 term, whose Fourier transform 
scales according to[| 



(O(x)O(O)) 



nI 



(5.6) 



9ymx 

This result passes the following important consistency test. The SYM in 2 dimensions 
with 16 supercharges have conformal fixed points in both UV and IR with central charges 
of order N 2 and N, respectively. Therefore, we expect the two point function of stress 
energy tensors to scale like N 2 /x 4 and N/x 4 in the deep UV and IR, respectively. Ac- 
cording to the analysis of [50|, we expect to deviate from these conformal behavior and 
cross over to a regime where supergravity calculation can be trusted. The cross over 
occurs at x = l/g Y AfVN and x = y/~N fgyM- At these points, the N scaling of ( |5.6|) and 



the conformal result match in the sense of the correspondence principle [[49 



5.2 Correlation functions from DLCQ 

The challenge then is to attempt to reproduce the scaling relation (|5.6|) , fix the numerical 
coefficient, and determine the detail of the cross-over behavior using SDLCQ. Ever since 
the original proposal [Q, the question of equivalence between quantizing on a light-cone 
and on a space-like slice have been discussed extensively. This question is especially 
critical whenever a massless particle or a zero-mode in the quantization is present. It is 
generally believed that the massless theories can be described on the light-cone as long 
as we take m — > as a limit. The issue of zero mode have been examined by many 
authors. Some recent accounts can be found in 



E6J, ML 10, B 



Generally speaking, 

supersymmetry seems to save SDLCQ from complicated zero-mode issues. We will not 
contribute much to these discussions. Instead, we will formulate the computation of the 
correlation function of stress energy tensor in naive DLCQ. To check that these results 
are sensible, we will first do the computation for the free fermions. Extension to SYM 



7 It is not difficult to show that for a generic p-brane, (O(x)O(O)) = N 5 -pg YA f~ 



50 



with 16 supercharges will be essentially straightforward, except for one caveat. In order 
to actually evaluate the correlation functions, we must resort to numerical analysis at the 
last stage of the computation. For the SYM with 16 supercharges, this problem grows 
too big too fast to be practical on desk top computer where the current calculations were 
performed. We can only provide an algorithm, which, when executed on an much more 
powerful computer, should reproduce (|5.6| ). Nonetheless, the fact that we can define a 
concrete algorithm seems to be a progress in the right direction. One potential pit-fall 
is the fact that the computation may not show any sign of convergence. If this is the 
case, or if it converges to a result at odds with Q5.6|) , we must go back and re-examine 
the issue of equivalence of forms and the issue of zero modes. 



The technique of DLCQ is reviewed by many authors |25|, |3(J so we will be brief here. 
The basic idea of light-cone quantization is to parameterize the space using light cone 
coordinates x + and x~ and to quantize the theory making x + play the role of time. In the 
discrete light cone approach, we require the momentum p_ = p + along the x~ direction 
to take on discrete values in units of p + /k where p + is the conserved total momentum 
of the system and k is an integer commonly referred to as the harmonic resolution. One 
can think of this discretization as a consequence of compactifying the x~ coordinate on a 
circle with a period 2L = 2irk/p + . The advantage of discretizing the light cone is the fact 
that the dimension of the Hilbert space becomes finite. Therefore, the Hamiltonian is a 
finite dimensional matrix and its dynamics can be solved explicitly. In SDLCQ one makes 
the DLCQ approximation to the supercharges and these discrete representations satisfy 
the supersymmetry algebra. Therefore SDLCQ enjoys the improved renormalization 
properties of supersymmetric theories. Of course, to recover the continuum result, we 
must send k to infinity and as luck would have it, we find that SDLCQ usually converges 
faster than the naive DLCQ. Of course, in the process, the size of the matrices will grow, 
making the computation harder and harder. 

Let us now return to the problem at hand. We would like to compute a general 
expression of the form 

F(x~,x + ) = (O(x~,x + )O(0,0)} . (5.7) 

In DLCQ, where we fix the total momentum in the x~ direction, it is more natural to 
compute its Fourier transform 

F(P„,x+) = -L(0(P_,x+)0(-P_,O)> . (5.8) 

This can naturally be expressed in a spectrally decomposed form 

F(P.,x + ) = J2^ Z (0\O(P-)\z}e-^ + {t \ O (-P^0)\0) ■ (5.9) 
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5.3 Correlator for Free Dirac Fermions 

Let us first consider evaluating this expression for the stress-energy tensor in the theory 
of free Dirac fermions as a simple example. The Lagrangian for this theory is 

C = iW$W-mti!V (5.10) 

where for concreteness, we take 7 = c 2 ,^ 1 = ia 1 and we take \1/ = 2~ 1//4 (^). In terms 
of the spinor components, the Lagrangian takes the form 

C = iijj*d+i/j + ix*d-x ~ "/s&V ~ ^**) • (5.11) 

Since we treat x + as time and since \ does not have any derivatives with respect to x + in 
the Lagrangian, it can be eliminated from the equation of motion, leaving a Lagrangian 
which depends only on ip: 

m 1 

£ = i< l p*d + <ip + i — VjriJ ■ (5-12) 

We can therefore express the canonical momentum and energy as 

P- = I dx~ iip*d^ip 

P+ = Jdx~ -!|V^V>- (5-13) 

In DLCQ, we compactify x~ to have period 2L. We can then expand ip and ip* in modes 

■0 = -j= {b{n)e~ s ^ x ~ + c/(-n)e^ ZLx ~] 

i/j* = -L (6(-n)e^*~ + (l(n)e-T 1 ") . (5.14) 



2L 

Operators b(n) and d{n) with positive and negative n are interpreted as a destruction 
and creation operators, respectively. In a theory with only fermions, it is customary to 
take anti-periodic boundary condition in order to avoid zero-mode issues. Therefore, n 
will take on half-integer values^. They satisfy the anticommutation relation 

{b(n), b(-m)} = {d(n), d(-m)} = 5„, m . (5.15) 

Now we are ready to evaluate ( |5.9|) in DLCQ. As a simple and convenient choice, we take 



0(-k) = - / </./■" (icO-c ■ ■i(0-c-)c)(- i T L - r . (5.16) 



3 In SDLCQ one must use periodic boundary condition for all the fields to preserve the supersymmetry. 
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which is the Fourier transform of the local expression for P with the total derivative 
contribution adjusted to make this operator Hermitian. Therefore, this should be thought 
of as the T ++ component of the stress energy tensor. For reasons that will become clear 
as we go on, this turns out to be one of the simplest things to compute. When acted on 
the vacuum, this operator creates a state 



r-(-*)|o> = I | 



n) b{-k + n)d{-ri)\0) . 



(5.17) 



Since the fermions in this theory are free, the plane wave states 



\n) = b(-k + n)d(-n)\0) 



(5.1* 



constitute an eigenstate. The spectrum can easily be determined by commuting these 
operators: 

I k k \ 

Ml\n) = 2P_P + \n) = m 2 - + In) (5.19) 

\n k — n J 

which is simply the discretized version of the spectrum of a two body continuum. All 
that we have to do now is calculate eigenstates of the actual theory we are interested in 
and to assemble these pieces into ( |5.9|) , but we can do a little more to make the result 
more presentable. The point is that since ( |5.9| ) is expressed in mixed momentum/position 
space notation in Minkowski space, the answer is inherently a complex quantity that is 
cumbersome to display. For the computation of two point function, however, we can go 
to position space by Fourier transforming with respect to the L variable. After Fourier 
transforming, it is straight forward to Euclideanize and display the two point function 
as a purely real function without loosing any information. To see how this works, let us 
write (|5.9f) in the form 



F(P_,x^ 



-(0\T ++ (k)\n) 



2LL 2 



e 2(hTL )x + 



(5.20) 



The quantity inside the absolute value sign is designed to be independent of L. Now, to 
recover the position space form of the correlation function, we inverse Fourier transform 
with respect to P = kir/L. 



Fix ,x 



T 2 

-(0\T ++ {k)\n) 

TV 



^(f) 1 



71 



kir \ L 



2(B3L) 



2n 2L L 2 
The integral over L can be done explicitly and gives 



(5.21) 



F(x , x + ) 



T 2 

-(0\T ++ {k)\n) 

TV 




Mi 



57T 2 fc 3 



K 4 (M n \f2 



x^x 



(5.22) 
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where K±(x) is the 4-th modified Bessel's function. We can now continue to Euclidean 
space by taking r 2 = 2x + x~ to be real and considering the quantity 




L 



7T 



(0\T ++ (k)\n) 



2 Mi 



•>7T 



A- 3 



K 4 {M n r) . (5.23) 



This is a fundamental result which we will refer to a number of times in this paper. It 
has explicit dependence on the harmonic resolution parameter k, but all dependence on 
unphysical quantities such as the size of the circle in the x~ direction and the momentum 



along that direction have been canceled. For the free fermion model, ( |5.23|) evaluates to 

iV M* {k-2n^ 




x+) -TV3^—k^ K ' iM " r) (5 ' 24) 



with M 2 given by (|5.19|) . The large k limit can be gotten by replacing n — > kx and 



-r J2 n ~~ *■ Jo dx- We recover the identical result using Feynman rules. For r <C m 1 , this 
behaves like 

JV 3(A;-2n) 2 N 




' ! ' = k ? 2vr2Pr4 "" 27T 2 r 4 ' (5 ' 25) 

5.4 Correlator for Supersymmetric Yang-Mills Theory with 16 
Supercharges. 

Finally, let us turn to the problem of computing the two point function of the T ++ op- 
erator for the SYM with 16 supercharges. Adopting light-cone coordinates and choosing 
the light-cone gauge will eliminate the gauge boson and half of the fermion degrees of 
freedom. The most significant change comes from the fact that the fields in this the- 
ory are in the adjoint rather than the fundamental representations and the theory is 
supersymmetric. This does not cause any fundamental problem in the DLCQ formu- 
lation of these theories. Indeed, the SDLCQ formulation of this || as well as many 
other related models with adjoint fields have been studied in the literature. The main 
difficulty comes from the fact that in supersymmetric theories low mass states such as 
tr[6(— n\)b{— ri2)b{— k+nx+ri2)]\0) with an arbitrary number of excited quanta, or "bits," 
appear in the spectrum. This means that for a given harmonic resolution k, the dimen- 
sion of the Hilbert space grows like exp(vfc), which is roughly the number of ways to 
partition k into sums of integers. 

The fact that the size of the problem grows very fast is somewhat discouraging from a 
numerical perspective. Nevertheless, it is interesting to note that DLCQ provides a well 
defined algorithm for computing a physical quantity like the two point function of T ++ 
that can be compared with the prediction from supergravity. In the following, we will 
show that this can be computed for the SYM theory by a straight forward application of 

(EM)- 
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As we found in section p] the momentum operator P + is given by 



P H 



dx tr 



(d-Xi) 2 + iu a d-U a 



(5.26) 



The local Hermitian form of this operator is given by 



T ++ (x) = tr 



1 



(d-X 1 ) 2 + - (iu a d^u a - i(d_u a )u a ) 



1... 



a = l...\ 



(5-27) 

where X and u are the physical adjoint scalars and fermions respectively, following the 
notation of ||. When discretized, these operators have the mode expansion 



X 



i -j 



<i 



^—[AS j (n)e-* m *-' L + jq i (-n)e 

n=l V n 

i oo 

= E [Bt j {n)e- i1 ' nx ~ /L + B^-ny™*' ~ /L 

4-^ n =l 



J7rn.a; /L 



In terms of these mode operators, we find 



T++(-A;)|0> 
fe-i 

E 



7T 

2L 



ra=l 



-y/n(k - n)Aij(-k + n)Aji(-n) + - - n\ S y (-/e + n)Bji(-n) 



(5.28) 

(5.29) 
|0). 



Therefore, (L/7r)(0|T ++ (— fe)|n) is independent of L and can be substituted directly into 
( ^.23| ) to give an explicit expression for the two point function. 

We see immediately that ( |5.29| ) has the correct small r behavior, for in that limit, 
( |5.29| ) asymptotes to (assuming n& = n/) 



x 



F(x ,x~ 



(5.30) 



N 2 v (3(k - 2nfn f 3n(k - n)n b \ _ N 2 (2n b + n f ) 
~Y V V 4vr 2 A; 2 r 4 + n 2 k 2 r 4 J ~ 4vr 2 r 4 



1 



which is what we expect for the theory of n& free bosons and n/ free fermions in the 
large k limit. 

Computing this quantity beyond the small r asymptotics, however, represents a 
formidable technical challenge. In we constructed the mass matrix explicitly and 
compute the spectrum for k = 2, k = 3, and k = 4. Even for these modest values of 
the harmonic resolution, the dimension of the Hilbert space was as big as 256, 1632, and 
29056 respectively (the symmetries of the theory can be used to reduce the size of the 
calculation somewhat). In figure |3], we display the results with the currently available 
values of k, except for the fact that we display the correlation function multiplied by a 
factor of 4Tf 2 r A /N 2 (2nb + n./), so that it now asymptotes to 1 (or in the logarithmic 
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Figure 3: (a) The spectrum as a function of 1/k and (b) the Log-Log plot of the correlation 
function (T ++ (x)T ++ (0)) (^+) N^t-inl+nA v - s - r * n ^ e unr ts where g\ M Nji{ = 1 for 



k = 3 and k = 4. 



Af 2 (2n 6 +n / ) 



scale) in the k —>■ oo limit. In this way any deviation from the asymptotic behavior 
1/r 4 is made more transparent. Note that with the values of the harmonic resolution k 
obtained at present, the spectrum in figure [|.a is far from resembling a dense continuum 
near M — 0. Clearly, we must probe much higher values of k before we can sensibly 
compare our field theory results with the prediction from supergravity. 



5.5 Supersymmetric Yang-Mills Theories with Less Than 16 
Supercharges 

The computation of the correlator for the stress energy tensor in the (8,8) model is limited 
by our inability to carry out the computation for large enough harmonic resolution. It is 
the (8,8) model which we are ultimately interested in solving in order to compare against 
the prediction of Maldacena's conjecture in the supergravity limit. Nevertheless, the 
computation of the correlation function can just as well be applied to models with less 
supersymmetry. We will conclude by reporting the results of such a computation. 

First, let us consider the theory with supercharges (1,1). This theory is argued not to 
exhibit dynamical supersymmetry breaking in |60], 0. We can also provide a physicist's 
proof that supersymmetry is not spontaneously broken for this theory by adopting the 
argument of Witten for the 2+1 dimensional SYM with Chern-Simons interaction [[84 



In [HI, the index of 2+1 dimensional SYM with gauge group SU(N) and 2 supercharges 
on R x T 2 was computed and was found to be non-vanishing for Chern-Simons coupling 
k 3 > N/2. If the period L of one of the circles in T 2 is sufficiently small, this theory 
is approximately the 2-dimensional SYM with (1,1) supersymmetry with gauge coupling 
! — g\j ' L and BF coupling k 2 
keeping g 2 and k 3 fixed. 



92 

L 



kuL [22 



Imagine approaching this theory by taking 
In this limit, k 2 — > in the units of g 2 so the limiting 
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Figure 4: (a) The spectrum as a function of 1/k and (b) the Log-Log plot of the correlation 

function (T ++ (x)T ++ (0)) ( |+ J j^-a^n +n ) V-S- r * n ^ ne un its where g\ M N /it = 1 for 
fc = 4...10. 



theory is that of pure SYM with (1,1) supersymmetry and a vanishing BF coupling. 
Choosing different values of k% corresponds to a different choice in the path of approach 
to this limit. If we chose k$ > N/2, we are guaranteed to have a non-zero index for 
finite L. This means that there will be a state with zero mass in the L — > limit also, 
indicating that supersymmetry is not spontaneously broken in this limit. On the other 
hand, the index is not a well defined quantity in the L — > limit, as a different choice of 
&3 will lead to a different value of the index in the L — > limit. In fact, the index can be 
made arbitrarily large by taking k% to be also arbitrarily large. This suggests that there 
are infinitely many states forming a continuum near m = 0. The index is therefore an ill 
defined quantity, akin to counting the number of exactly zero energy states on a periodic 
box as one takes the volume to infinity. 



This theory is also believed not to be confining [15[ [16| and is therefore expected to 
exhibit non-trivial infra-red dynamics. 

The SDLCQ of the 1+1 dimensional model with (1,1) supersymmetry was solved in 
p4] , [|, and we apply these results directly in order to compute ( p.23| ) . For simplicity, we 
work to leading order in the large N expansion. The spectrum of this theory for various 
values of k, and the subsequent computation of (|5.23|) is illustrated in figure f|.a. 

The spectrum of this theory at finite k, illustrated in figure f|.a, consists of 2k — 2 
exactly massless states^, accompanied by large numbers of massive states separated by 
a gap. The gap appears to be closing in the limit of large k however. We have tried 
extrapolating the mass of the lightest massive state as a function of 1/k by performing a 
least square fit to a line and a parabola, giving the extrapolated value of M 2 T[ /g\ M N = 
1.7 and M 2 -k / g\ M N = —0.6, suggesting indeed that at large k, the gap is closed. This 
is consistent with the expectation that the spectrum is that of a continuum starting at 



i.e. k — 1 massless bosons, and their superpartners. 
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Figure 5: ■^J2 n \^(^\T ++ (k)\n)\ 2 v.s. k from states with M\n) = 0. This quantity 
determines the coefficient of the 1/r 4 asymptotic tail of the correlation function in the 
large r limit for the (1,1) model. 

M = discussed earlier, although one must be careful when the order of large N and 
large k limits are exchanged. At finite N, we expect the degeneracy of 2k — 2 exactly 
massless states to be broken, giving rise to precisely a continuum of states starting at 
M = as expected. 

In the computation of the correlation function illustrated in figure [|.b, we find a 
curious feature that it asymptotes to the inverse power law c/r 4 for large r. This be- 
havior comes about due to the coupling (0|T ++ |n) with exactly massless states \n). The 
contribution to Q5.231 ) from strictly massless states are given by 
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(a) The spectrum as a function of 1/k and (b) the Log-Log plot of the correlation 

(T++(.x)T++(0)> (^) 2 ^g;+n / ) v -s- r in the units where & M N/n = 1 for 
6. 



x 



X" 1 



F(x , x + ) 



IX 



'■(0\T ++ (k)\n) 



M: 



8ir 2 k 3 



T^K A (M n r) 



(5.31) 



M„=0 



r 2 

-<0|T ++ (£;)|n> 

7T 



6 



j\y 



n =o /c 3 7r 2 r 4 



We have computed this quantity as a function of 1/k and extrapolated to 1/k — > by 
fitting a line and a parabola to the computed values for finite 1/k. The result of this 
extrapolation is illustrated in figure [|. The data currently available suggests that the 
non-zero contribution from these massless states persists in the large k limit. 

Let us now turn to the model with (2,2) supersymmetry. The SDLCQ version of this 
model was solved in fill . The result of this computation can be applied to (|5.23| ). The 
result is summarized in figure |6|. This model appears to exhibit the onset of a gapless 
continuum of states more rapidly than the (1,1) model as the harmonic resolution k is 
increased. Just as we found in the (1,1) model, this theory contains exactly massless 
states in the spectrum. These massless states appear to couple to T ++ |0) only for k 
even, and the overlap appears to be decreasing as k is increased. We believe that this 
model is likely to exhibit a power law behavior c/r 1 for 7 > 4 for the T ++ correlator for 
r ^> gYM\N in the large N limit. Unfortunately, the existing numerical data do not 
permit the reliable computation of the exponent 7. 
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6 Bound States of Three Dimensional Supersym- 
metric Theory. 

Recently, there has been considerable progress in understanding the properties of strongly 



coupled gauge theories with supersymmetry |75|, |74)| . In particular, there are a number 
of supersymmetric gauge theories that are believed to be inter-connected through a web 
of strong-weak coupling dualities. Although these dualities provide a deep insight into 
the dynamics of gauge theory at strong and weak couplings, they do not usually give 
much information about the spectrum of bound states at intermediate values of coupling 
constant g. The prominent exception is so called BPS states whose mass is protected by 
supersymmetry and thus stays the same for all values of g. An interesting new possibility 
for analytical treatment of bound state problem in SYM^ is based on the duality between 
SYM and supergravity proposed by Maldacena |nj as discussed in the previous section. 
This idea was also exploited in [^] to get the glueball spectrum of three dimensional 
theory and the results agree with lattice calculations. 

However it would be interesting to solve the bound state problem for SYM theory 
directly, starting from the first principles of quantum field theory. As we have seen in 
the previous sections the solution can be found for the various two dimensional theories 
by means of applying discrete light cone quantization. Evidently, it would be desirable 
to extend these DLCQ/SDLCQ algorithms to solve higher dimensional theories. One 
important difference between two dimensional and higher dimensional theories is the 
phase diagram induced by variations in the gauge coupling. The spectrum of a 1 + 1 
dimensional gauge theory scales trivially with respect to the gauge coupling, while a 
theory in higher dimensions has the potential of exhibiting a complex phase structure, 
which may include a strong-weak coupling duality. It is therefore interesting to study 
the phase diagram of gauge theories in D > 3 dimensions. 

Towards this end, we consider three dimensional SU(iV) M = 1 super- Yang-Mills 
compactified on the space-time R x S 1 x S 1 . In particular, we compactify the light- 
cone coordinate x~ on a light-like circle via DLCQ, and wrap the remaining transverse 
coordinate x x on a spatial circle. By retaining only the first few excited modes in the 
transverse direction, we are able to solve for bound state wave functions and masses 
numerically by diagonalizing the discretized light-cone supercharge. We show that the 
supersymmetric formulation of the DLCQ procedure - which was studied in the context 
of two dimensional theories extends naturally in 2 + 1 dimensions, resulting in an exactly 
supersymmetric spectrum. 

6.1 Light-Cone Quantization and SDLCQ 

We wish to study the bound states of M = 1 super- Yang-Mills in 2 + 1 dimensions. Any 
numerical approach necessarily involves introducing a momentum lattice - i.e. parton 
momenta can only take on discretized values. The usual space-time lattice explicitly 
breaks supersymmetry, so if we wish to discretize our theory and preserve supersymmetry, 
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then a judicious choice of lattice is required. 

In 1 + 1 dimensions, it is well known that the light-cone momentum lattice induced 
by the DLCQ procedure preserves supersymmetry if the supercharge rather than the 
Hamiltonian is discretized |64|, [7|. In 2 + 1 dimensions, a supersymmetric prescription is 
also possible. We begin by introducing light-cone coordinates x ± = (x° ± x 1 )/^, and 
compactifying the x~ coordinate on a light-like circle. In this way, the conjugate light- 
cone momentum k + is discretized. To discretize the remaining (transverse) momentum 
k L = k 2 , we may compactify x x = x 2 on a spatial circle. Of course, there is a significant 
difference between the discretized light-cone momenta k + , and discretized transverse 
momenta k±; namely, the light-cone momentum k + is always positive 10 , while k± may 
take on positive or negative values. The positivity of k + is a key property that is exploited 
in DLCQ calculations; for any given light-cone compactification, there are only a finite 
number of choices for k + - the total number depending on how finely we discretize the 
momenta^]. In the context of two dimensional theories, this implies a finite number of 
Fock states ||68|| . 

In the case we are interested in - in which there is an additional transverse dimension - 
the number of Fock states is no longer finite, since there are an arbitrarily large number of 
transverse momentum modes defined on the transverse spatial circle. Thus, an additional 
truncation of the transverse momentum modes is required to render the total number of 
Fock states finite, and the problem numerically tractable^. In this work, we choose the 
simplest truncation procedure beyond retaining the zero mode; namely, only partons with 
transverse momentum k± = 0, ±^ will be allowed, where L is the size of the transverse 
circle. 

Let us now apply these ideas in the context of a specific super- Yang-Mills theory. We 
start with 2 + 1 dimensional M = 1 super- Yang- Mills theory defined on a space-time with 
one transverse dimension compactified on a circle: 

S = Jd 2 xJ dx ± tr(~F^F^ + M^D^). (6.1) 

After introducing the light-cone coordinates x ± = 4= (x° ± x 1 ), decomposing the spinor 
\1/ in terms of chiral projections - 



„5 1 „,5 



*=-$-*> X=^I73T* ( 6 - 2 ) 



and choosing the light-cone gauge A + = 0, the action becomes 



S = I dx + dx / <ir_i_tr 



i(cLA-) 2 + D + 0«9_0 + itpD + ij+ 



10 Since we wish to consider the decompactificd limit in the end, we omit zero modes. This is a 
necessary technical constraint in numerical calculations. 

n The 'resolution' of the discretization is usually characterized by a positive integer K, which is called 
the 'harmonic resolution' [p8L p3|; for a given choice of K, the light-cone momenta k + are restricted to 
positive integer multiples of P + /K, where P + is the total light-cone momentum of a state 

12 This truncation procedure, which is characterized by some integer upper bound, is analogous to the 
truncation of k + imposed by the 'harmonic resolution' K . 
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+ixd-X + -j^D L (t) + -= 



Jl> 



(6.3) 



A simplification of the light-cone gauge is that the non-dynamical fields A and \ may 
be explicitly solved from their Euler-Lagrange equations of motion: 



A- 



X 



9 T= 9 
1 



iJ=i(^,M + 2#) 



(6.4) 



y/2d- 



D ± ip. 



These expressions may be used to express any operator in terms of the physical degrees 
of freedom only. In particular, the light-cone energy, P~ , and momentum operators, 
P + ,P ± , corresponding to translation invariance in each of the coordinates x and x± 
may be calculated explicitly: 



P' 
P- 



dx / dx±tr 
Jo 

dx~ / dx±tr 
Jo 



\d_(f)) 2 + iipd^ip] , 


(6.5) 


.-^-l^. 


(6.6) 


d-4>d±(J) + hpd±ip] . 


(6.7) 



The light-cone supercharge in this theory is a two component Majorana spinor, and may 
be conveniently decomposed in terms of its chiral projections: 



Q- 



2 1/4 / dx- / dx ± tr [0<9_^ - ^d. 



2 3/4 / dx~ / dxxtx 



2d ± <M>+g(i[<f>,d-<i>]+2iH;)—il> 



(6.8) 
(6.9) 



The action ( |6.3| ) gives the following canonical (anti) commutation relations for propagat- 
ing fields at equal x + \ 

4>ij(x~,x±), d-(/)ki(y~, y±) = -ti(x~ - y~)5(x± - y±) (Su8j k - —8ijS k ij ,(6.10) 
{^ij(x~,x ± ),ip k i(y~,y ± )} = -5{x~ -y~)6(x± - y±) \5 a 5 jk - —5ij5ki\ ■ (6.11) 
Using these relations one can check the supersymmetry algebra: 
{Q + ,Q+} = 2v / 2> + , {Q-,Q-} = 2V2P-, {Q + ,Q~} = -4P±. (6.12) 

We will consider only states which have vanishing transverse momentum, which is 
possible since the total transverse momentum operator is kinematicaQ On such states, 

13 Strictly speaking, on a transverse cylinder, there are separate sectors with total transverse momenta 
2im/L; we consider only one of them, n = 0. 
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the light-cone supercharges Q + and Q~ anti-commute with each other, and the super- 
symmetry algebra is equivalent to the M = (1,1) supersymmetry of the dimensionally 
reduced (i.e. two dimensional) theory p3j| . Moreover, in the P± = sector, the mass 
squared operator M 2 is given by M 2 = 2P + P~ . 

As we mentioned earlier, in order to render the bound state equations numerically 
tractable, the transverse momentum of partons must be truncated. First, we introduce 
the Fourier expansion for the fields <fi and ip, where the transverse space-time coordinate 
x x is periodically identified: 



<f>ij(Q,X ,x ± 

1 

V2^L 
ipij(0,x~,x± 



E 



00 dk A 



n-^=—oo 



V2k+ 



dij(k + , n L )e 1 



k+ 



X — 1 1 X 



lx n t (U+ vMJk+x +ii 



+ a!j i (k + ,n )e 



-x± 



v n J -=— 00 



b ij (k + ,n x )e- ik+x -- i '^^ + b} i (k + ,n ± )e lk+ *- +i2 - 2 ^ 



Substituting these into the (anti) commutators (|6.11|) , one finds 



t ( n + 



a*j(p ,nx),a\ k (q + ,m± 



5(p + - q + )5 n±jm± ( SuSjk - jrSijSik 



{bij{p + , n±), b\ k (q + , mi)} = 5(p + - q + )5 n± , m± \8u8 jk - —5^5^ 
The supercharges now take the following form: 

g+ ==i2 i/4 J2 )T M difev^[6l i (A;,n- L )o ii (A;,n- L )-oI i (A;,n- L )6y(A;,n- L )' 

Q- 

4 



(6.13) 
(6.14) 

(6.15) 






y J dk—= a\j(k,n ± )bij(k,n x ) — bl j (k,n ± )a i j(k,n ± ) 
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i2 _1 / 4 (7 roc 
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(6.16) 
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We now perform the truncation procedure; namely, in all sums over the transverse mo- 
menta nj- appearing in the above expressions for the supercharges, we restrict summation 
to the following allowed momentum modes: nj- = 0, ±1. More generally, the truncation 
procedure may be defined by \nf\ < N max , where N max is some positive integer. In this 
work, we consider the simple case N max = 1. Note that this prescription is symmetric, 
in the sense that there are as many positive modes as there are negative ones. In this 
way we retain parity symmetry in the transverse direction. 

How does such a truncation affect the supersymmetry properties of the theory? Note 
first that an operator relation [A, B] = C in the full theory is not expected to hold in 
the truncated formulation. However, if A is quadratic in terms of fields (or in terms of 
creation and annihilation operators), one can show that the relation [A, B] = C implies 



[Afri B tr \ — C, 



tr 



for the truncated operators A tr ,B tr , and C tr - In our case, Q + is quadratic, and so the 
relations {Qtr,Qt} — ^V^Ptr an d {QtnQtr} — are true in the P± — sector of the 
truncated theory. The {Qtr, Qtr} however is not equal to 2\[2P^.. So the diagonalization 
of {Qfr, Qtr} will yield a different bound state spectrum than the one obtained after 
diagonalizing 2\f2P tT .. Of course the two spectra should agree in the limit N max — > oo. 
At any finite truncation, however, we have the liberty to diagonalize any one of these 
operators. This choice specifies our regularization scheme. 

Choosing to diagonalize the light-cone supercharge, however, has an important ad- 
vantage: the spectrum is exactly super symmetric for any truncation. In contrast, the 
spectrum of the Hamiltonian becomes supersymmetric only in the N max —>■ oo limitf^. 

To summarize, we have introduced a truncation procedure that facilitates a numer- 
ical study of the bound state problem, and preserves supersymmetry. The interesting 
property of the light-cone supercharge Q~ [Eqn fl6.16"D l is the presence of a gauge coupling 



constant as an independent variable, which does not appear in the study of two dimen- 
sional theories. In the next subsection, we will study how variations in this coupling 
affects the bound states in the theory. 

6.2 Numerical Results: Bound State Solutions 

In order to implement the DLCQ formulation of the bound state problem - which is 
tantamount to imposing periodic boundary conditions x~ = x~ + 2ttR []63] - we simply 



restrict the light-cone momentum variables ki appearing in the expressions for Q + and 
Q~ to the following discretized set of momenta: l-^P + ,j^P + ,j^P + ,..., i. Here, P + 
denotes the total light-cone momentum of a state, and may be thought of as a fixed 
constant, since it is easy to form a Fock basis that is already diagonal with respect to the 



operator P + [B8|. The integer K is called the 'harmonic resolution', and 1/K measures 



14 If one chooses anti-periodic boundary conditions in the x coordinate for fermions, then there is no 
choice; one can only diagonalize the light-cone Hamiltonian. See [E9| for more details on this approach. 
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the coarseness of our discretization. The continuum limit is then recovered by taking 
the limit K — > oo. Physically, 1/K represents the smallest positive unit of longitudinal 
momentum-fraction allowed for each parton in a Fock state. 
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Figure 7: Plot of bound state mass squared M 2 in units 16ir 2 N/L 2 as a function of the 
dimensionless coupling < g' < 2, defined by (g 1 ) 2 = g 2 NL/16ir 3 , at N = 1000 and K = 5. 
Boson and fermion masses are identical. 

Of course, as soon as we implement the DLCQ procedure, which is specified unam- 
biguously by the harmonic resolution K, and cut-off transverse momentum modes via 
the constraint \nf\ < N max , the integrals appearing in the definitions for Q + and Q~ 
are replaced by finite sums, and so the eigen-equation 2P + P~|\f') = M 2 |\I>) is reduced 
to a finite matrix diagonalization problem. In this last step we use the fact that P~ is 
proportional to the square of the light-cone supercharge^ Q ■ in the present work, we 
are able to perform numerical diagonalizations for K = 2,3,4 and 5 with the help of 
Mathematica and a desktop PC. In Figure |7|, we plot the bound state mass squared M 2 , 
in units 16tt 2 N/L 2 , as a function of the dimensionless coupling g' = g\/NL/4:ir 3 ^ 2 , in the 
range < g' < 2. We consider the specific case N = 1000, although our algorithm can 
calculate masses for any choice of N, since it enters our calculations as an algebraic vari- 
able. Since there is an exact boson-fermion mass degeneracy, one needs only diagonalize 
the mass matrix M 2 corresponding to the bosons. For K = 5, there are precisely 600 
bosons and 600 fermions in the truncated light-cone Fock space, so the mass matrix that 

15 Strictly speaking, P~ = -j=(Q~") 2 is an identity in the continuum theory, and a definition in the 
compactified theory, corresponding to the SDLCQ prescription p\ M . 
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Figure 8: Plot of bound state mass squared M 2 in units 16tt 2 N/L 2 as a function of the 
dimensionless coupling < g' < 10, defined by {g') 2 = g 2 NL/lQir 3 , at N = 1000 and K = 5. 
Note the appearance of a new massless state at strong coupling. 

needs to be diagonalized has dimensions 600 x 600. At K = 4, there are 92 bosons and 
92 fermions, while at K = 3, one finds 16 bosons and 16 fermions. 

In Figure H, we plot the bound state spectrum in the range < g' < 10. It is apparent 
now that new massless states appear in the strong coupling limit g' — > oo. 

An interesting property of the spectrum is the presence of exactly massless states 
that persist for all values of the coupling g'. For K — 5, there are 16 such states (8 
bosons and 8 fermions). At K = 4, one finds 8 states (4 bosons and 4 fermions) that are 
exactly massless for any coupling, while for K = 3, there are 4 states (two bosons and two 
fermions) with this property. We will have more to say regarding these states in the next 
subsection, but here we note that the structure of these states become 'string-like' in the 
strong coupling limit. This is illustrated in Figure |9], where we plot the 'average length' (or 
average number of partons) of each of these massless states^. This quantity is obtained 
by counting the number of partons in each Fock state that comprises a massless bound 
state, appropriately weighted by the modulus of the wave function squared. Clearly, at 
strong coupling, the average number of partons saturates the maximum possible value 
allowed by the resolution - in this case 5 partons. The same behavior is observed at 
lower resolutions. Thus, in the continuum limit K — > oo, we expect the massless states 



16 The 'noisiness' in this plot for larger values of g' reflects the ambiguity of choosing a basis for the 
eigen-space, due to the exact mass degeneracy of the massless states. 
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Figure 9: Plot of average length for the eight massless bosonic states as a function of the 
dimensionless coupling g' , defined by (g') 2 = g 2 NL/16ir 3 , at N = 1000 and K = 5. Note that 
the states attain the maximum possible length allowed by the resolution K = 5 in the limit of 
strong coupling. 

in this theory to become string-like at strong coupling. 

One interesting property of the model studied here is the manifest Af = (1, 1) super- 
symmetry in the P x = momentum sector, by virtue of the supersymmetry relations 
( |6.12| ). Moreover, if we consider retaining only the zero mode nf~ = 0, then the light-cone 
supercharge Q~ for the 2 + 1 model is identical to the 1 + 1 dimensional Af = (1, 1) su- 
persymmetric Yang-Mills theory studied in |64], |5], |J, after a rescaling by the factor 1/ g' . 
(This is equivalent to expressing the mass squared M 2 in units g 2 N/7r, where g = g/yL. 
The quantity g is then identified as the gauge coupling in the 1 + 1 theory) We may 
therefore think of the additional transverse degrees of freedom in the 2 + 1 model, repre- 
sented by the modes n 1 = ±1, as a modification of the 1 + 1 model. A natural question 
that follows from this viewpoint is: How well does the 1 + 1 spectrum approximate the 
2 + 1 spectrum after performing this rescaling? Before discussing the numerical results 
summarized in Table |2|, let us first attempt to predict what will happen at small coupling 
g'. In this case, the coefficients of terms in the rescaled Hamiltonian P~ that correspond 
to summing the transverse momentum squared l^ -1 ! 2 of partons in a state will be large. 
So the low energy sector will be dominated by states with n 1 - = 0. i.e. those states that 
appear in the Fock space of the Af = (1,1) model in 1 + 1 dimensions. This is indeed 
supported by the results in Table 0. 

For large coupling g', however, it is clear that the approximation breaks down. In 
fact, one can show that the tabulated masses in the rescaled 2 + 1 model tend to zero in 
the strong coupling limit, which eliminates any scope for making comparisons between 
the two and three dimensional models. 
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Comparison Between 1 + 1 and 2 + 1 Spectra 


- 


1 + 1 Model 


Rescaled 2 + 1 Model 


K 


- 


g> = .01 


g' = -i 


g< = 1.0 


K = 5 


15.63 

18.23 
21.8 


15.5 
17.6 
21.3 


15.17 
17.9 
21.7 


3.7 
3.5 
3.2 


K = 4 


18.0 
21.3 


17.99 
21.3 


17.6 
21.0 


3.56 
3.1 


K = 3 


20.2 


20.2 


19.8 


3.1 



Table 2: Values for the mass squared M 2 , in units g 2 N/Ti, with g 2 = g 2 /L, for bound states in 
the dimensionally reduced N = (1, 1) model, and the 2 + 1 model studied here. The quantity g 
is identified as the gauge coupling in the 1 + 1 model. We set K = 3, 4 and 5, and N = 1000. 
Note that the comparison of masses between the 1 + 1 model, and the (re-scaled) 2 + 1 model 
is good only at weak coupling g' '. 



Thus, the non-perturbative problem of solving dimensionally reduced models in 1 + 1 
dimensions can only provide information about bound state masses in the corresponding 
weakly coupled higher dimensional theory. 

6.3 Analytical Results: The Massless Sector 

In the previous subsection we presented the results of studying the bound state problem 
using numerical methods. In performing such a study we conveniently chose the simplest 
nontrivial truncation of the transverse momentum modes; namely, n± = 0, ±1. Surpris- 
ingly, such a simple scheme provided many interesting insights concerning the massless 
and massive sector. In particular we see that there are three types of massless states; 
those that are massless only at g = or g = oo (but not both), and those that are massless 
for any value of the coupling. In this subsection, we will analyze only the massless sector 
of the theory, and show that the observed properties of the spectrum with the truncation 
n 1 - = 0, ±1 also persists if we include higher modes: n 1 - = 0, ±1, ±2, . . . , ±N max . We 
therefore consider the model with supercharges given by fl6.15| ) and ( |6.16| ), and restrict 
summation of transverse momentum modes via the constraint {n^l < N max . 

For states carrying positive light-cone momentum, P + is never zero, and so massless 
states must satisfy the equation -P~|4 7 ) = 0, which, using the relation P = -^(Q - ) 2 , 
and hermiticity of Q~ , reduces to 

Q~|tf) = 0. (6.17) 



This is the equation we wish to study in detail. 
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We begin with an analysis of the weak coupling limit of the theory. This limit means 
that the dimensionless coupling constant is small: i.e. gyL <C 1. We will consider the 
strong-weak coupling behavior of the theory on a cylinder with fixed circumference L so 
it is convenient to choose the units in which L = 1 for this discussion. The supercharge 
( |6.16| ) consists of two parts: one is proportional to the coupling and the other is coupling- 
independent: 

Q~ = Q± + gQ- (6.18) 

So at g = 0, the equation ( |6.171 ) reduces to Q±\*&) = 0, which means that |\&) may be 



viewed as a state in the Fock space of the two dimensional M = (1, 1) super Yang-Mills 
theory, which may be obtained by dimensional reduction of the 2 + 1 theory. Thus the 
massless states at g = are states with any combination of a/(k,0) and tf(k,0) modes, 
and no partons with nonzero transverse momentum. 

What happens with these massless states when one switches on the coupling? To 
answer this question, we need some information about the behavior of states as functions 
of the coupling. We assume that wave functions are analytic in terms of g at least in the 
vicinity of g — 0. This means that in this region any massless state |\1/) may be written 
in the form: 

oo 

|*> = £<?», (6.19) 

n=0 



where states \n) are coupling independent. Then using relation ( |6.18| ), the ^-dependent 



equation ( 6.171 ) may be written as an infinite system of relations between different \n): 



Q±\0) = 0, (6.20) 

Q±\n)+Q\n-1) = Q, n>0. (6.21) 

The first of these equations was already used to exclude partons carrying non-zero trans- 
verse momentum, which is a property of the massless bound states at zero coupling. The 
second equation is non-trivial. Let us consider two different subspaces in the theory. The 
first of these subspaces consists of states with no creation operators for transverse modes 
which we will label 1. The other is the complement of this space in which the operator 
Q± is invertible and we label this space 2. Equation ( |6.21| ) defines the recurrence relation 
when Q \ n — 1) is in subspace 2: 

\n) = -Ql 1 (Q\n-l)\ 2 ), (6.22) 

The consistency condition is that projection of Q\n — 1) in subspace 1 is zero, 

Q|n-l}| =0. (6.23) 

This condition implies that not all states of the two dimension theory, g = , may be 
extended to such states at arbitrary g using ( |6.22| ). Taking n = 1, ( |6.23| ) implies that 
|0) is a massless state of the dimensionally reduced theory. The numerical solutions, of 
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course, show this correspondence between the 2 + 1 and 1 + 1 |>4], |], massless bound 
states. Starting from a massless state of the two dimensional theory, and we construct 
states \n) using ( |6.22j ), and for which ( |6.23|) is always satisfied. Then \fy) may be found 
from summing a geometric series: 

OO -| 

l*> = E(-^: 1 ^)"|0) = — ^- 10). (6.24) 

n=0 1 + flV_|_ V 

So, starting from the massless state of the two dimensional M = (1,1) model, one can 
always construct unique massless states in the three dimensional theory at least in the 
vicinity of g — 0. 

The state (|6.24|) turns out to be massless for any value of the coupling: 



Q-|*> = Q±(l + gQl'Q)- ^TTdO) = Q±\0) = 0, (6.25) 

1 + gQ ± Q 

though the state itself is dependent on g. Thus, we have shown that massless states 
of the three dimensional theory, at nonzero coupling, can be constructed from massless 
states of the corresponding model in two dimensions. All other states containing only 
two dimensional modes can also be extended to the eigenstates of the full theory. But 
such eigenstates are massless only at zero coupling. Assuming analyticity, one can then 
show that their masses grow linearly at g in the vicinity of zero. Such behavior also 
agrees with our numerical results. 

To illustrate the general construction explained above we consider one simple example. 
Working in DLCQ at resolution K = 3 we choose a special two dimensional massless 
statesman: 

|0) =tr(o t (l,0)o t (2,0))|uoc). (6.26) 

Then in the SU(iV) theory we find: 

Q\0) = ^[tr(ot(l,0)(6t(l ) -l)ot(l,l)-ot(l,l)6t(l ) -l)+ 

+ 6 f (l, 1)^(1, -1) - o f (l, -l)6 f (l, 1)))] \vac), (6.27) 

|1) = -g 1 1 g|0>=- i ^^=(at(l,0)at(l,-l)at(l,l)- 

- a + (l, 0)a f (l, l)a f (l, — 1)) \vac) (6.28) 

Q\l) = 0. (6.29) 

The last equation provides the consistency condition ( |6.23|) for n = 2, and it also shows 
that for this special example we have only two states |0) and |1), instead of a general 
infinite set. The matrix form of the operator 1 + gQ^Q in the |0), |1) basis is 

i + gQl 1 Q=( 1 \ 9 ) = (l i) ■ (6-30) 



17 The state |0) denotes a massless state, while \vac) represents the light-cone vacuum. 
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Then the solution of ( |6.24| ) is 

|tf) = |0) +g\l) = tr(o f (l, 0)a f (2, 0))\vac) + (6.31) 

+ ££27f (° t(1 ' ° )at(1 ' 1)at(1 ' _1) ~ ° t(1 ' ° )at(1 ' ~ 1)at(1 ' 1} ) lVaC) - 

This state was observed numerically, and the dependence of the wave function on the 
coupling constant is precisely the one given by the last formula. 

In principle, we can determine the wave functions of all massless states using this 
formalism. Our procedure has an important advantage over a direct diagonalization of the 
three dimensional supercharge. Firstly, in order to find two dimensional massless states, 
one needs to diagonalize the corresponding supercharge |63]]. However, the dimension of 



the relevant Fock space is much less than the three dimensional theory (at large resolution 
K, the ratio of these dimensions is of order (N max + l) aK , a ~ 1/4). The extension of 
the two dimensional massless solution into a massless solution of the three dimensional 
theory requires diagonalizing a matrix which has a smaller dimension than the original 
problem in three dimensions. Thus, if one is only interested in the massless sector of the 
three dimensional theory, the most efficient way to proceed in DLCQ calculations is to 
solve the two dimensional theory, and then to upgrade the massless states to massless 
solutions in three dimensions. 

Finally, we will make some comments on bound states at very strong coupling. Of 
course, we have states (|6.24|) which are massless at any coupling, but our numerical 



calculation show there are additional states which become massless at g = 00 (see Figure 
H). To discuss these state it is convenient to consider 

Q- = -Q± + Q (6.32) 

9 

instead of Q~, and perform the strong coupling expansion. Since we are interested 
only in massless states, the absolute normalization doesn't matter. We repeat all the 
arguments used in the weak coupling case: first, we introduce the space 1* where Q can 
not be inverted, and its orthogonal complement 2*. Then any state from 1* is massless 
at g = 00, but assuming the expansion 

00 1 
!*> = £-»* (6-33) 

n=o y 



at large enough g, one finds the analogs of ( |6.22| ) and ( |6.23|) : 



\n)* = -Q- 1 {Q±\n-l)*\ 2 .), (6.34) 

Q±\n-l)*\ 1 . = 0. (6.35) 

As in the small coupling case, there are two possibilities: either one can construct all 
states \n)* satisfying the consistency conditions, or at least one of these conditions fails. 
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The former case corresponds to the massless state in the vicinity of g — oo, which can 
be extended to the massless states at all couplings. The states constructed in this way 



- and ones given by ( 6.24 ) - define the same subspace. In the latter case, the state is 
massless at g = oo, but it acquires a mass at finite coupling. There is a big difference, 
however, between the weak and strong coupling cases. While the kernel of Q± consists 
of "two dimensional" states, the description of the states annihilated by Q is a nontrivial 
dynamical problem. Since the massless states can be constructed starting from either 
g = or g = oo, we don't have to solve this problem to build them. If, however, one 
wishes to show that massless states become long in the strong coupling limit (there is 
numerical evidence for such behavior - see Figure ^), the structure of 1* space becomes 
important, and we leave this question for future investigation. 

Conclusion. 

In these lectures we have reviewed some of the progress in the application of discrete 
light cone quantization to supersymmetric systems. Studying such systems is especially 
interesting because the cancellation between bosonic and fermionic loops make these 
theories much easier to renormalize than the models without supersymmetry. Although 
we didn't need this advantage when considering two dimensional systems, it becomes 
crucial in higher dimensions. From this point of view it is desirable to have exact SUSY 
in discretized theories to simplify the renormalization in DLCQ. 

while we are still far from the point of solving the bound state problem in three and 
four dimensional theories, we can already make some statements about these theories. 
For example in section § we described the vacuum structure of SYM 2 +i on a cylinder. 
The reason for this is that only zero modes contribute to such structure, thus studying 
the theory dimensionally reduced to 1 + 1 provide all the necessary information. As we 
saw in section g, two dimensional models can also be used to determine the behavior of 
bound states at weak coupling in three dimensions and to count the exact massless states. 
We performed such counting only for (1,1) theory, the case of (2,2) supersymmetry |~1| 



and the even more interesting case of the (8,8) theory || which is known to have a mass 
gap have not been addressed. 

Let us now mention a few of the immediate challenges for SDLCQ. First of all is 
the extention of the numerical results of section ^ to higher resolution and thus to test 
the Maldacena's conjecture. The only problem here is the limits in one's computing 
resources and the the speed of the algorithems one uses. Impovements here will also 
help use to extend our analysis of higher dimensional system and to larger values of 
transverse truncation. Unfortunately the transverse truncation that we have achieved 
so far does not provide much information about behavior of the spectrum as a function 
of transverse resolution and was used mainly for illustration of the general concepts. 
Another consideration is that our approach studies theories on a light like cylinder and 
thus our theories may have as aspects in which they are different than theories in infinite 
spacetime or on a space like cylinder. 
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